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ABSTRACT 

Accurate techniques that allow for the derivation of the spatial density in star formation regions are 
rare. A technique that has found application for the derivation of spatial densities in Galactic star 
formation regions utilizes the density-sensitive properties of the K-doublet transitions of formaldehyde 
(II2CO). In this paper, we present an extension of our survey of the formaldehyde lip — In (A — 6.2 cm) 
and 2ii — 2i2 (A = 2 .1cm) K-doublet transitions of II2CO in a sample of 56 starburst systems 



(Mangum et al. 20081. We have extended the number of galaxies in which both transitions have 
been detected from 5 to 13. We have improved our spatial density measurements by incorporating 
kinetic temperatures based u pon NH3 measurem ents of 11 of the galaxies with a total of 14 velocity 
components in our sample (Man gum et al.|2013 ). Our spatial density measurements lie in a relatively 
narrow range from lO-*-^ to lO'^^'^ cm This implies that the Schmidt-Kennicutt relation between 
L/fl and yidense'- (1) is an indication of the dense gas mass reservoir available to form stars, and 
(2) is not directly dependent upon a higher average density driving the star formation process in the 
most luminous starburst galaxies. We have also used our H2CO measurements to derive two separate 
measures of the dense gas mass which are generally smaller, in many cases by a factor of 10^ — 10^, 
than those derived using HCN. This disparity suggests that II2CO traces a denser, more compact, 
component of the giant molecular clouds in our starburst galaxy sample. 

We also report measurements of the rotationally-excited A = 6.3cm ^ni/2^ — 1/2 state of OH and 
the Hllla radio recombination line taken concurrently with our H2CO lio — In measurements. 
Subject headings: ISM: molecules, galaxies: ISM, galaxies: nuclei, galaxies: starbursts, radio lines: 
galaxies 



1. INTRODUCTION 

Studies of the abundant CO molecule toward external 
galaxies have revealed the large scale structure of the 
molecular mass in these objects. Unfortunately, high op- 
tical depth tends to limit the utility of the abundant CO 
molecule as a probe of the spatial density and kinetic 
temperature in dense gas environments, quantities nec- 
essary t o assess the possib ility and conditions of star for- 
mation (jLiszt et al.j2010|). Emission from less-abundant, 
higher-dipole moment molecules is better-suited to the 
task of deriving the spatial density and kinetic tem- 
perature within dense gas environments. A variety of 
molecules which trace dense molecular environments in 
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mainly nearby galaxies have been stud ied 



see the re- 
for early 



views of H enkel, Baan, fc Mauersberger ( 1991 
work) and |Omont| (|2007| . 

Formaldehyde (H2CO) has proven to be a reliable den- 
sity and kinetic temperatur e probe in Galactic molecu- 



lar clouds ( Mangum & Wootten 1993 
19931 pnsburg et al.||2Dllff 
the HzCO lio 



Mangum et al. 

Existing measurements ot 
2i2 K-doublet transitions 



111 and 2ii 

in a w ide variety of gala,xies b y Baan et aT] ([r986 , 199(5 



1993|, and |Araya et al.| ( |2004| have mainly cohcentratec 

111 



on measurements ot the lio— In transition. By collecting 
a consistent set of measurements of both K-doublet tran- 
sitions we are able to derive the spatial density within 
the extragalactic star formation regions which comprise 
this study. Using the unique density selectivity of the K- 
doublet transitions of H2CO we have measured the spa- 
tial density in a sample of galaxies exhibiting starburst 
phenomena and/or high infrared luminosity (Mangum 
et al.|[2008 ). In the present work we have nearly tripled 
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the number of galaxies within which the spatial density 
can be derived from five to thirteen. In eleven of these 
galaxies with a total of fourteen velocity components we 
have further improved the quality of our spatial density 
determination by applying kinetic t emperatures derived 
from NII3 measurements (Mangum et al. 2013). In [2^ 
we discuss the specific properties of the H-jcU molecule 
that make it a good probe of spatial density. Jpjpresents a 
summary of our observations; Qour H2CO, OH, Hllla, 
and continuum emission measurement results; f|5] analy- 
ses of our H2CO, OH, and Hllla measurements, includ- 
ing Large Velocity Gradient (LVG) model fits and dense 
gas mass calculations based on our H2CO measurements. 

2. FORMALDEHYDE AS A HIGH DENSITY PROBE 

The ubiquitous and readily-observed formaldehyde 
(H2CO) molecule is a reliable probe of the high den- 
sity environs of star formation regions. Pervasive in the 
interstellar medium, H2CO exhibits a moderate abun- 
dance which appears to vary by less than an order- 
of-magni tude within star form ation regions in our own 



Galaxy (Mangum et al. [1990 ). As a slightly asymmet- 
ric rotor molecule, most H2CO rotational transitions are 
split into two sub levels. These energy levels are desig- 
nated by a total angular momentum quantum number, 
J, the projection of J along the symmetry axis for a lim- 
iting prolate symmetric top, K_i, and the projection of J 
along the symmetry axis for a limiting oblate symmetric 
top, K-|_i. This splitting leads to two types of transi- 
tions: the high-frequency A J ~ 1, AK_i = 0, AK+i 
= — 1 "P-branch" transitions and the lower- frequency 
AJ = 0, AK_i = 0, AK+i = ±1 "Q-branch" transi- 
tions, popularly known as the "K-doublet" transitions 



the lower 
- 5i5 transi- 



see discuss ion in Mangum & Wootten 1993 McGauley 
et al. 2011). The P-branch transitions are only seen in 
emission in regions where n(iJ2) ^ 10'' cm^^. In the K 
doublet transitions, for n{H2) < 10^'^ cm 
energy states of the lio — In through 5i4 
tions b ecome overpqpulated due to a collisional s election 
effect ( [Evans et"aL]|l975[ [Garrison et al.[[1975l ). This 
overpopulation cools the J< 5 K-doublets to excitation 
temperatures lower than that of the cosmic microwave 
background, causing them to appear in absorption. For 
n(i72) ^ 10^'^ cm~'^ and an assumed kinetic temperature 
of 40 K, this collisional pump is quenched and the J < 
5 K-doublets are driven into emission over a wide range 
of kinetic temperatu res and abundances (see Figure 1 in 



Mangum et al. 



2008). 



Measurements of the spatial density in galactic and 
extragalactic star formation regions often rely upon only 
order-of-magnitude or critical density estimates. The col- 
lisional excitation of H2CO provides a direct translation 
to density that is fairly kinetic temperature independent. 
The use of H2CO as a density probe offers unprecedented 
improvement to these more approximate density deter- 
minations, allowing for an actual measurement of this 
important physical quantity. 

3. OBSERVATIONS 

The measurements reported here were made using the 
National Radio Astronomy Observatory (NRACj^ Green 

^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 



Bank Telescop e (GBT) during the p eriods 2006/06/07- 
14 (reported in [Mangum et al.j2008[ ), 2007/01/16-04/22, 
2007/11/05-11/26, 2008/04/11-09/22, and 2011/05/04- 
07/16. Using exactly the same observi ng setup as was 
used for the observations presented in [Mangum et al.| 
(i2008|), single-pointing measurements w ere obtain ed of 
theTio - 111 (4.829660 GHz; [T ucker et ^[l97l| ) and 
2ii -2i2 (14.488479 GHz; [Tucker et al., 1 9 72p K-doublet 



transitions of H2CO, the Hllla radio recombination 
line (RRL) at 4.744183 GHz, and two (i^ = 1 - 
and 1 — 1) of the three 6 cm lines of the rotationally- 
excited ^Hi/a J = 1/2 state of OH at 4.7506 56(3) and 
4.765562(3) GHz ( Radfordj [19681 |Lovas[|19"92| toward a 
sample of 56 infrared-luminous and/or starburst galax- 
ies (Table [1]). Our galaxy sample was chosen to represent 
both galaxies with measured H2GO em ission or absorp- 
tion ( |Baan et"aL 1993 Araya et al. 2004 ) and other galax- 
ies with substantial molecular emission deduced from 



HCN measurements ( Gao & Solomon 2004a ) . The single 



beam H2CO lio - In, Hllla, and OH ^Hi/2J = 1/2 
transitions were measured simultaneously using 4 spec- 
tral windows each with 50 MHz of bandwidth sampled 
by 16384 channels. All of the 4.8 GHz measurements 
{9b = 153") utilized the position switching technique 
with reference position located 30 arcmin west in az- 
imuth from each galaxy position. Galaxy and reference 
integrations were 3 minutes each at 4.8 GHz. The H2CO 
2ii — 2i2 transition was measured using a dual-beam 
{9b — 51"; beam separation 330" in cross-elevation) 
receiver over 50 MHz of bandwidth sampled by 16384 
channels. The dual-beam system at 14.5 GHz allowed for 
both position switching or beam nodding. Galaxy and 
reference integrations were 2 minutes each at 14.5 GHz. 
The correlator configurations produced a spectral chan- 
nel width of 3.052 kHz, which is approximately 0.2 and 
0.08 km s^^ at 4.8 and 14.5 GHz, respectively. 

For eighteen of the galaxies listed in Table [T] mea- 
surements of one or both K-doublet H 2CO transitions 
were presented in Mangum et al. (2008). Additional ob- 
servations and/or updated analysis (i.e. modified spec- 
tral smoothing) of the K-doublet H2CO measurements 
of these eighteen galaxies are included in the results pre- 
sented here. 

To calibrate the intensity scale of our measurements, 
several corrections need to be considered: 

Opacity:: Historical opacity estimates based on atmo- 
spheric model calculations using ambient pressure, 
temperature, and relative humidity measurements 
indicated that r at 4.8 and 14.5 GHz was ~ 0.01 
and 0.05 during our observations (assuming ele- 
vation > 30 degrees). The opacity corrections 
exp(ro csc(i5L)) are < 1.02 and < 1.05, respec- 
tively. 

Flux:: Assuming point-source emission, one can 
use the current relation (derived from point- 
source radiometric continuum measurements; 
[2008 



Maddalena 



for the aperture efficiency 
rjA = 0.71 exp (0.0163i^(GHz))^^ to convert an- 
tenna temperature to flux density. At 4.8 and 14.5 

ment by Associated Universities, Inc. 
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TABLE 1 

EXTRAGALACTIC FORMALDEHYDE SURVEY SOURCE LiST 



Galaxy^ 


a 


6 






d 


Glassification° 




(J2000) 


(J2000) 


(km s^i) 


(Mnc) 


(K) 




NGC55 


00:14:54.5 


-39:11:19 


129 


1.5 + 0.1 




SB(s)m: edge-on 


NOG 253 


00:47:33.1 


-25:17:18 


251 


3.4 + 0.2 


34 


SAB(s)c 


IC 1623 


01:07:47.2 


-17:30:25 


6028 


80.7 + 5.7 


39 


LIRG 


NGG520 


01:24:35.3 


+03:47:37 


2281 


30.5 + 2.1 


38 


Pec, Pair, Sbrst 


NGG598 


01:33:54.0 


+30:40:07 


-179 


0.9 + 0.3 




SA(s)cd 


NGG660 


01:43:01.7 


+13:38:36 


856 


12.2 + 0.9 


37 


SB(s)a:pec, HII LINER 


IR01418+1651 


01:44:30.5 


+17:06:09 


8101 


109.7+7.7 




LIRG 


NGG695 


01:51:14.9 


+22:34:57 


9769 


130.2 + 9.1 


34 


SO pec:LIRG 


Mrkl027 


02:14:05.6 


+05:10:24 


9061 


120.8 + 8.5 


37 


LLIRG 


NGG891 


02:22:33.4 


+42:20:57 


529 


9.4 + 0.7 


28 


SA(s)b?:sp 


NGG925 


02:27:16.9 


+33:34:35 


553 


9.3 + 0.7 




SAB(s)d 


NGG1022 


02:38:32.7 


-06:40:39 


1503 


19.2+ 1.3 


39 


SB(s)a 


NGG1055 


02:41:45.2 


+00:26:35 


996 


13.4 + 0.9 


29 


SBb:sp:Sy2 LINER 


Maffei 2 


02:41:55.1 


+59:36:15 


-17 


3.1 ± 0.2 




SAB(rs)bc 


NGG 1068 


02:42:40.7 


-00:00:48 


1136 


15.2+ 1.1 


40 


SA(rs)b:Syl/2 


UGG 02369 


02:54:01.8 


+14:58:15 


9262 


124.8 + 8.8 




DBL 


NGG 1144 


02:55:12.2 


-00:11:01 


8750 


115.3 + 8.1 


32 


S pec 


NGG 1365 


03:33:36.4 


-36:08:25 


1652 


21.5+ 1.5 


32 


SBb(s)b 


IR 03359+1523 


03:38:47.1 


+15:32:53 


10507 


142.1 ± 10.0 




LIRG 


IG 342 


03:46:49.7 


+68:05:45 


31 


3.8 + 0.3 


30 


SAB(rs)cd 


NGG 1614 


04:33:59.8 


-08:34:44 


4847 


64.2 + 4.5 


46 


SB(s)c:pec 


VIIZw31 


05:16:46.4 


+79:40:13 


16290 


220.8+ 15.5 


34 


NGG 2146 


06:18:37.7 


+78:21:25 


918 


16.7+ 1.2 


38 


SB(s)ab:pec 


NGG 2623 


08:38:24.1 


+25:45:17 


5535 


79.4 ± 5.6 




LIRG 


Arp55 


09:15:55.1 


+44:19:55 


11957 


164.3+ 11.6 


36 


Pair 


NGG 2903 


09:32:10.1 


+21:30:02 


556 


7.4 + 0.5 


29 


SB(s)d 


UGG 5101 


09:35:51.6 


+61:21:11 


11810 


164.3+ 11.5 


36 


Syl.5, LINER 


M82 


09:55:52.2 


+69:40:47 


203 


5.9 + 0.4 


45 


IO,Sbrst 


M82SW 


09:55:50.0 


+69:40:43 


203 








NGG 3079 


10:01:57.8 


+55:40:47 


1150 


20.7+ 1.5 


32 


SB(s)c, LINER 


IR 10173+0828 


10:19:59.9 


+08:13:34 


14716 


206.7+ 14.5 




Sbrst 


NGG 3227 


10:23:30.7 


+19:52:46 


1111 


20.3+ 1.4 




SAB(s):pec 


NGG 3627 


11:20:15.0 


+12:59:30 


727 


6.5 + 0.5 


30 


SAB(s)b:LINER 


NGG 3628 


11:20:17.2 


+13:35:20 


847 


8.5 + 0.6 


30 


Sb:pec:sp 


NGG 3690 


11:28:32.2 


+58:33:44 


3121 


48.5 + 3.4 




Merger 


NGG 4631 


12:42:08.0 


+32:32:29 


606 


7.6 + 0.5 


30 


SB(s)d 


NGG 4736 


12:50:53.0 


+41:07:14 


308 


4.8 + 0.3 




SA(r)ab;Sy2;LINER 


Mrk231 


12:56:14.2 


+56:52:25 


12642 


178.1 ± 12.5 




SA(rs)c:pec 


NGG 5005 


13:10:56.2 


+37:03:33 


946 


19.3+ 1.4 


28 


SAB(rs)bc 


IG 860 


13:15:04.1 


+24:37:01 


3866 


53.8 + 3.8 




Sa, Sbrst 
SA(s)bc:pec 


NGG 5194 


13:29:52.7 


+47:11:43 


463 


9.1 ± 0.6 




M83 


13:37:00.9 


-29:51:57 


518 


4.0 + 0.3 


31 


SAB(s)c 


Mrk273 


13:44:42.1 


+55:53:13 


11324 


160.5+ 11.2 


48 


LINER 


NGG 5457 


14:03:12.6 


+54-20-57 


241 


6.2 ± 0.4 




SAB(rs)cd 


IR 15107+0724 


15:13:13.1 


+07:13:27 


3897 


61.9 + 4.4 




Sbrst 


Arp 220 


15:34:57.1 


+23:30:11 


5434 


82.9 + 5.8 


44 


Pair, Sbrst 


NGG 6240 


16:52:59.0 


+02:24:02 


7339 


108.8 + 7.6 


41 


IO:pec, LINER, Sy2 


IR 17208-0014 


17:23:21.9 


-00:17:00 


12834 


183.0 + 12.8 


46 


ULIRG 


IR 17468+1320 


17:49:06.7 


+13:19:54 


4881 


74.1 ± 5.2 




LIRG 


NGG 6701 


18:43:12.4 


+60:39:12 


3950 


59.1 ± 4.1 


32 


SB(s)a 


NGG 6921 


20:28:28.8 


+25:43:24 


4399 


63.1 ± 4.4 


34 


SA(r)0/a 


NGG 6946 


20:34:52.3 


+60:09:14 


48 


5.5 + 0.4 


30 


SAB(rs)cd 


IG5179 


22:16:09.1 


-36:50:37 


3447 


48.8 + 3.4 


33 


SA(rs)bc 


NGG 7331 


22:37:04.1 


+34:24:56 


821 


14.4+ 1.0 


28 


SA(s)b;LINER 


NGG 7479 


23:04:56.6 


+12:19:22 


2385 


33.7+2.4 


36 


SB(s)c;LINER 


IR 23365+3604 


23:39:01.3 


+36:21:09 


19330 


262.5+ 18.4 


45 


ULIRG 


Mrk331 


23:51:26.7 


+20:35:10 


5422 


74.9 ± 5.2 


41 


LIRG 



IR = IRAS (Infrared Astronomical Satellite) throughout this paper. 
Heliocentric velocity drawn from the literature. 

NED^ Hubble flow distance corrected for Virgo cluster, Great Attractor, and Shapley supercluster. For NGG 598 no Hubble 
flow di stance available, so NED "redshift-independent" distance assumed. 

Prom [Gao fc Solomon] |2004b| , who used IRAS 60 and 100/im dust continuum emission ratios with an assumed dust emissivity 
oc v-P with /3 = 1.5. 
" From NED, Sbrst = starburst galaxy. 

^ The NASA/IPAG Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, Galifornia Institute of Tech- 
nology, under contract with the National Aeronautics and Space Administration. 
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GHz this yields tja — 0.71 and 0.67, respectively. 
For elevation 90 degrees and zero atmospheric 
opacity, T^/S = 2.846//^ = 2.0 and 1.9 for 4.8 
and 14.5 GHz, respectively. For elevation > 30°, 
Ta/S = 1.97 and 1.80 at 4.8 and 14.5 GHz. These 
are the K/Jy calibration factors used to convert 
our spectra to flux density assuming point-source 
emission. Measurements of the standard flux 
calibration sources 3C 48 and 3C 286 yielded 
Ta/S = 1.95 ± 0.04 and 1.83 ± 0.04 at 4.8 and 
14.5 GHz, consistent with the standard empirical 
values. Note also that since the opacity correction 
is small, = TAexp(Aro) ~ ^ ~ T^, where 
rii = 0.99 for the GET. Using r?Jh ~ 1.327?^, we 
can write the main beam brightness temperature 



Galaxy Structure:: Several of our galaxies are sus- 
pected to have high-density structure (measured 
with high-dipolc moment molecules like HCN) on 
scales approaching the size of our beam at 14.5 GHz 
(51"). The beam coupling correction necessary to 
account for structure in our 14.5 GHz measure- 
ments, relative to our point-source assumption, is 
given by: 



fcoupl 



92 



ing 



92 



(1) 



This correction factor is less than 20% for 9s < 23". 
With the exception of M82, none of the galaxies 
in our sample have measured dense molecular gas 
structure larger than ^ 20". We therefore assume 
that, with the exception of M8^ aU of the H2CO 
emission reported in these measurements is from 
structures smaller than the primary beams of our 
measurements . 

Absolute Amplitude Calibration:: The GET abso- 
lute amplitude calibration is reported to be accu- 
rate to 10-15% at all frequencies, limited mainly 
by temporal drifts in the noise diodes used as ab- 
solute amplitude calibration standards. As noted 
above, measurements of the standard flux calibra- 
tion sources 3C48 and 3C286 yielded T^/S = 
1.95±0.04 and 1.83±0.04 at 4.8 and 14.5 GHz, sug- 
gesting that our absolute flux calibration is good 
to ~ 5%. Relative calibration between our 4.8 and 
14.5 GHz measurements is assumed to be ~ 5%, 
which produces a negligible impact on density mea- 
surements obtained from H2CO line ratios. 

4. RESULTS 
4.1. H2CO 

Table[2] lists our H2CO lio — In and 2ii — 2i2 mea- 
surement results, including updated results for the eigh- 
te en galaxies in our pres ent sample that were presented 



and/or updated analysis of these eighteen galaxies are 
included in the analysis presented in this work, we list 
these updated H2CO measurement results in italics in 
Table [2j H2CO measurements for thirty-eight galaxies 
and one galaxy offset position (M 82SW) are completely 
new measurements. The new H2CO measurements pre- 
sented in this paper have nearly tripled the number of 
galaxies within which both H2CO K-doublet transitions 
have been detected from five to thirteen. For each detec- 
tion we list the peak intensity, heliocentric central line 
velocity, velocity width (FW(ZI/HM)), and integrated 
intensity derived from direct channel-by-channel integra- 
tion of and gaussian fits to each of the line profiles. Un- 
certainties expressed as one-sigma are listed in parenthe- 
ses for each quantity in Tables [2j 

Spectra for the galaxies detected in only the lio — In or 
2ii — 2i2 transition are displayed in Figure [T] NGC 253, 
NGC660, Maffei2, IC342 (ah in Figure NGC 2146, 
M82, M82SW, NGC 3079 (all in Figure 3 1 NGC 3628, 
IC860, M83, IR 15107-^0724 (aU in Figured, Arp220, 
and NGC 6946 (both in Figure ^ were all aetected in 
both H2CO transitions. These spectra have been gaus- 
sian smoothed to the spectral channel widths quoted 
m Table [2] to both increase the channel-specific signal- 
to-noise ratio of our measurements and more closely 
match other mole cular spectral line mea surements of 
Gao"fc Solomon|[2004a ). 



these galaxies (i.e. 

4.2. OH and Hllla 

Table [3] lists and Figures |6] and [7] show our measured 
OH^Hi/aJ 1/2F = 1-0 and 1-1 [E^ ~ 182K above 
ground) and Hllla radio recombination line results. We 
report new detections of these OH and Hllla transitions 
toward nine galaxies, while updated analysis is presented 
for OH and Hllla detections quo ted for four galaxies 



presented in [Mangum et al. ( 2008 ) . For each OH and 



Hllla detection we list the same quantities in Table |3] 
as derived for our H2CO spectra in Table [2j 

Toward most of our galaxies we did not detect the OH 
or Hllla transitions. The RMS noise values for our OH 
4765 and 4750 MHz and Hllla measurements of the 
undetected galaxies are listed in Table |4j 

Rotationally excited OH has been used as a tracer of 
the molecu l ar environment within AGN (|Henkel et al. 



1986 



[19871 [Henkel fc Wilson[ [T990{ [Impellizzeri et al. 



2006[ ). AGN come in two main types: those with (type 
1) and those without (type 2) broad optical atomic line 
emission. In the unified scheme of active galactic nuclei, 
all AGN are intrinsically similar, with the observed dif- 
ferences in type due to orientation/observer perspective. 
In the framework of a common paradigm a significant 
column density of molecular material, in the form of a 
parsec-scale torus, obscures our view of the AGN in type 
2 objects. Attempts to detect the obscuring molecular 
material using measurements of molecular abso rption or 



emission h ave yielded few c onfirmati ons (e.g. Schmelz 
et al 1986} [Baa.n fc Haschick 1992; St aveley-Smith et al. 
19921). While the existence of molecular AGN tori is 



Mangum et al.[ ( |2008[ ). Since additional observations '^^H ffrmed by H2O megamaser emission in some objects 



^ As was discussed iii |Mangum et ah] ( 2008| l, the source coupling 
correction for the extended nature ot M 82 results in a 20% decrease 
in the measured H2CO lio ~ lll/2ii — 2i2 K-doublet transition 
ratio, which leads to a corresponding increase in the derived spatial 
density. 



( [Lo[[2005[ Reid et al. 2009), other tracers like OH m 
the ground rotational transitions at 1.7 GHz tra ce gas 
at larger galactocentric distance (e.g. [Pih lstrb m et al. 
2001| |Klockner et al.[[2003{ [Klockner fc Baan||2004] ) Be- 



fore concluding that parsec-scale tori are rarely molec 
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TABLE 2 
H2CO Measurements* 



Galaxy 


Transition*^ 


Fit" 






FW(HM/ZI)'=''' 
(km s-i) 


/T::idv« 








(mK) 


(km s"^) 


(mK km s-^) 



NGC55 


2ii — 2i2 


DIO 


(0.8) 








NGC253 


iio ~ ^11 


D15 


-32.8(1.8) 


230.1 


349.1 


-5773.7(299.5) 






G15 


-34.3(1.3) 


227.0(3.0) 


163.0(7.0) 


-5947.9(339.4) 




2ii — 2i2 


D15 


-20.2(0.4) 


237.8 


242.9 


-2542.3(51.6) 






G15 


-11.4(1.0) 


169.1(3.1) 


58.3(6.3) 


-708.2(99.5) 






G15 


-18.7(0.6) 


252.6(2.9) 


93.1(6.3) 


-1847.8(136.4) 


IC 1623 


lio — 111 


DIO 


(1.3) 








NGC 520 


iio — ill 


D20 


-4.4(0.9) 


2301.2 


215.4 


-449.1(93.9) 






G20 


-4.0(0.4) 


2286.6(4.8) 


106.8(11.6) 


-451.9(64.3) 




2ii - 2i2 


DIO 


(0.5) 








NGC 598 


lio — 111 


DIO 


(0.7) 










2ii — 2i2 


DIO 


(0.3) 








NGC 660 


iici — ill 


D30 


-3.2(0.6) 


1041.6 


976.6 


-1511.9(261.7) 






G30 


-2.2(0.2) 


923.3(47.3) 


765.0(135.5) 


-1758.2(366.4) 




2n - 2i2 


D30 


-3.7(0.4) 


871.6 


1060.8 


-1280.8(223.5) 






G30 


— 2.2(0.2) 


754.3(19.1) 


547.2(47.8) 


— 1257.6(140.3) 


IR01418+1651 


lio — 111 


D20 


(0.6) 










2ii — 2i2 


DIO 


(0.3) 








NGC 695 


lio — 111 


DIO 


(0.8) 










2ii — 2i2 


DIO 


(0.5) 








Mrk 1027 


lio ~ 111 


DIO 


(1.3) 








NGC 891 


iio — ill 


D20 


-2.8(0.8) 


476.3 


434.8 


-538.9(160.1) 






G20 


-1.7(0.3) 


535.0(48.8) 


413.6(145.2) 


-726.6(280.4) 




2ii — 2i2 


DIO 


(0.3) 








NGC 925 


lio — 111 


DIO 


(1.4) 








NGC 1022 


lio ~ 111 


DIO 


(1.2) 










2ii — 2i2 


DIO 


(0.7) 








NGC 1055 


lio — 111 


DIO 


(1.2) 










2ii — 2i2 


DIO 


(0.3) 








Maffei 2 


lio — 111 


D20 


-5.5(0.9) 


-64.0 


354.5 


-756.5(144.7) 






G20 


-2.1(0.6) 


-98.1(45.5) 


127.5(85.2) 


-284.9(208.7) 






G20 


-4.5(1.0) 


24.7(17.1) 


102.1(30.3) 


-493.2(180.2) 




2ii — 2i2 


D20 


— 1.4(0.2) 


-15.5 


303.1 


— 220.5(24.7) 






G20 


— 0.9(0.2) 


-90.7(12.8) 


78.9(29.0) 


-78.1(33.2) 






G20 


— 1.4(0.2) 


18.3(9.1) 


91.3(21.8) 


— 136.0(36.9) 


NGC 1068 


lio — 111 


DIO 


(1.1) 










2ii — 2i2 


DIO 


(0.4) 








UGC 02369 


lio — 111 


DIO 


(1.1) 








NGC 1144 


lio — 111 


D20 


-6.9(0.7) 


8838.7 


1214.4 


-2251.6(417.6) 






G20 


— 2.1(0.7) 


8281.8(8.6) 


51.6(20.1) 


— 116.8(60.4) 






G20 


-6.2(0.4) 


8517.4(5.4) 


158.7(13.2) 


— 1044.8(113.0) 






G20 


-1.8(0.2) 


9088.9(39.0) 


609.0(119.6) 


-1186.0(275.9) 




2ii — 2i2 


DIO 


(0.4) 








NGC 1365 


lio — 111 


D20 


— 2.1(0.4) 


1637.6 


398.4 


—410.0(77.1) 






G20 


-1.6(0.3) 


1608.1(21.8) 


266.9(57.9) 


-449.3(121.8) 




2ii — 2i2 


DIO 


(2.0) 








IR 03359+1523 


lio — 111 


DIO 


(1.3) 








IC 342 


iio — ill 


DIO 


-5.6(1.7) 


22.1 


116.7 


-391.1(94.8) 






GIO 


-5.2(0.8) 


25.0(4.9) 


67.9(11.5) 


-372.2(83.4) 




2 11 - 2i2 


DIO 


-2.8(0.4) 


23.5 


109.8 


-138.8(22.7) 






GIO 


-2.6(0.2) 


27.4(2.1) 


51.1(4.9) 


-139.1(17.6) 


NGC 1614 


lio — 111 


DIO 


(1.5) 








VIIZw31 


lio — 111 


DIO 


(1.1) 










2ii — 2i2 


DIO 


(0.3) 








NGC 2146 


lio ~ 111 


D20 


-1.9(0.3) 


883.9 


436.0 


-483.2(68.7) 






G20 


-1.7(0.2) 


867.3(13.4) 


239.8(36.4) 


-430.3(77.3) 






G20 


-1.6(0.4) 


1038.5(5.3) 


49.6(14.5) 


-83.8(31.4) 




2ii — 2i2 


D20 


-1.7(0.1) 


869.4 


569.1 


-366.0(19.6) 






G20 


-1.2(0.1) 


824.3(8.9) 


231.8(23.3) 


-296.8(35.2) 






G20 


-0.6(0.1) 


1046.9(12.3) 


109.5(29.4) 


-65.3(22.1) 


NGC 2623 


lio — 111 


DIO 


(1.0) 










2ii — 2i2 


DIO 


(0.4) 








Arp55 


iio — ill 


DIO 


(2.6) 








NGC 2903 


iio ~ ill 


DIO 


(0.3) 










2ii — 2i2 


DIO 


(0.2) 








UGC 05101 


lio — ill 


D20 


(0.6) 










2ii — 2i2 


DIO 


(0.3) 








M82 


iio — ill 


D20 


-22.8(1.1) 


231.7 


453.8 


-4600.1(237.0) 






G20 


-14.7(0.9) 


124.6(4.1) 


98.1(9.0) 


-1531.6(167.6) 






G20 


-21.7(0.7) 


275.0(3.3) 


135.3(8.0) 


-3128.0(212.8) 




2ii - 2i2 


D20 


-6.0(0.4) 


234.6 


465.5 


-1440.9(92.7) 






G20 


-5.7(0.6) 


110.2(5.8) 


97.4(13.1) 


-595.2(99.4) 






G20 


-4.7(0.4) 


278.6(10.3) 


177.1(26.6) 


-883.5(151.2) 
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Galaxy 


Transition^ 


Fit" 


T\ 




FW(HM/ZI)°'' 










CmK') 


(km s^^) 


(km s"^) 


(mK km s~^) 


M82SW 


lio — 111 


D20 


-18.2(1.8) 


228.0 


424.4 


-4040.0(320.0) 






G20 


-11.2(1.4) 


132.0(12.6) 


116.4(24.6) 


-1384.6(342.9) 






G20 


-18.0(1.1) 


282.7(9.1) 


144.2(19.2) 


-2767.9(406.9) 




2ii — 2i2 


D20 


-7.7(0.4) 


225.2 


464.8 


-1632.9(98.1) 






G20 


-7.3(0.5) 


116.4(3.5) 


96.4(7.9) 


-743.3(80.2) 






G20 


-4.7(0.3) 


274.9(9.3) 


185.4(22.7) 


-925.2(125.0) 


NGC 3079 


iio — ill 


G05 


+ 1.0(0.4) 


1006.8(12.0) 


55.2(33.6) 


56.0(41.4) 






G05 


-2.0(0.2) 


1186.7(11.4) 


200.9(35.0) 


-436.1(87.9) 






G05 


-0.9(0.2) 


1501.9(25.6) 


213.5(78.8) 


-210.6(90.3) 




2ii — 2i2 


G05 


-5.0(0.6) 


1017.8(3.4) 


56.6(8.2) 


-300.9(57.0) 






G05 


-27.6(0.9) 


1115.8(0.4) 


25.5(1.0) 


-748.3(37.4) 


IR 10173+0828 


iio — ill 


DIO 


(2.2) 








NGC 3227 


lio ~ 111 


DIO 


(2.1) 








NGC 3627 


lio — 111 
2ii — 2i2 


DIO 
DIO 


(0.8) 
(0.3) 








NGC 3628 


iio — ill 


DIO 


-8.5(1.1) 


826.4 


384.0 


-1295.6(191.7) 






GIO 


-4.1(0.9) 


725.9(11.5) 


87.7(25.5) 


-381.7(137.8) 






GIO 


-5.8(0.6) 


875.4(11.4) 


153.8(29.0) 


-954.9(206.4) 




2ii — 2i2 


DIO 


-2.5(0.4) 


900.7 


210.8 


-209.4(43.2) 






GIO 


-1.7(0.2) 


886.5(6.4) 


119.1(16.0) 


-217.1(37.6) 


NGC 3690 


lio — 111 
2ii — 2i2 


DIO 
DIO 


(0.8) 
(0.4) 








NGC 4631 


2ii — 2i2 


DIO 


(0.3) 








NGC 4736 


2ii — 2i2 


DIO 


(0.8) 








Mrk231 


lio — 111 
2ii — 2i2 


DIO 
DIO 


(0.7) 
(0.4) 








NGC 5005 


2ii — 2i2 


DIO 


(0.7) 








IC 860 


iio — ill 


DIO 


+5.5(1.2) 


3921.1 


315.8 


501.5(185.3) 






GIO 


+3.6(0.6) 


3889.1(9.3) 


117.5(22.0) 


449.5(111.3) 




2ii — 2i2 


DIO 


-1.2(0.2) 


3879.8 


215.0 


-128.9(17.3) 






GIO 


-1.1(0.1) 


3877.9(3.6) 


108.8(8.6) 


-130.9(13.6) 


NGC 5194 


lio ~ 111 


D20 


-1.7(0.3) 


369.7 


434.5 


-222.8(67.3) 






G20 


-0.8(0.2) 


403.9(30.9) 


267.7(81.1) 


-237.2(89.8) 




2ii — 2i2 


DIO 


(0.3) 








M83 


iio — ill 


DIO 


-3.0(0.8) 


510.0 


181.0 


-167.6(67.1) 






GIO 


-1.7(0.5) 


522.6(15.6) 


104.3(39.2) 


-192.3(93.1) 




2ii - Zi2 


DIO 


-1.5(0.3) 


494.6 


152.1 


-117.6(21.1) 






GIO 


-1.1(0.2) 


507.2(8.7) 


109.2(24.0) 


-129.9(34.9) 


Mrk 273 


lio — 111 


DIO 


-4.7(1.5) 


11521.7 


171.6 


-216.6(120.3) 






GIO 


-2.5(0.8) 


11505.0(15.5) 


95.3(38.2) 


-253.0(131.2) 




2ii — 2i2 


DIO 


(0.4) 








NGC 5457 


lio ~ 111 
2ii — 2i2 


DIO 
DIO 


(0.9) 
(0.4) 








IR 15107+0724 


iiQ — ill 


DIO 


+5.3(1.0) 


3888.1 


391.9 


714.0(190.1) 






GIO 


+3.7(0.3) 


3888.4(8.0) 


182.0(19.1) 


713.5(98.5) 




2ii - 2i2 


DIO 


-3.7(0.4) 


3912.3 


368.6 


-506.7(75.5) 






GIO 


-2.7(0.2) 


3908.6(5.7) 


179.3(13.7) 


-508.9(51.2) 


Arp 220 


iio — ill 


D20 


+5.5(0.7) 


5483.5 


531.2 


1323.3(171.3) 






G20 


+5.1(0.4) 


5352.3(7.0) 


154.1(18.2) 


842.1(116.7) 






G20 


+3.2(1.0) 


5487.2(6.7) 


59.5(19.9) 


200.4(92.0) 






G20 


+ 1.7(0.4) 


5596.0(38.8) 


175.6(86.4) 


316.9(168.8) 




2ii - 2i2 


D20 


-4.1(0.4) 


5405.0 


544.5 


-1012.2(89.7) 






G20 


-3.5(0.9) 


5318.0(24.7) 


136.9(29.6) 


-505.4(167.7) 






G20 


-3.2(0.7) 


5443.1(25.5) 


137.6(42.7) 


-474.5(181.5) 






G20 


-0.5(0.2) 


5607.4(29.1) 


85.4(76.7) 


-41.1(42.9) 


NGC 6240 


iio ~ ill 


D20 


-1.9(0.6) 


7309.8 


392.8 


-354.3(101.8) 






G20 


-1.5(0.3) 


7309.5(26.7) 


260.4(70.2) 


-401.1(136.7) 




2ii — 2i2 


DIO 


(0.4) 








IR 17208-0014 


lio — 111 
2ii — 2i2 


DIO 
DIO 


(1.2) 
(0.9) 








liX 1- 1 ^OOt^ J-OZU 


^10 111 


DIO 


f 1 5) 








NGC 6701 


lio — 111 
2ii — 2i2 


DIO 
DIO 


(1.4) 
(0.5) 








NGC 6921 


lio — 111 


D20 


-3.4(0.6) 


3960.7 


609.9 


-1012.7(182.8) 






G20 


-2.8(0.3) 


3995.4(20.1) 


364.4(52.0) 


-1087.9(196.5) 




2ii — 2i2 


DIO 


(0.4) 








NGC 6946 


iio — ill 


DIO 


-2.7(0.7) 


51.2 


275.9 


-269.2(95.5) 






GIO 


-1.7(0.5) 


-11.6(6.8) 


44.4(16.4) 


-78.5(37.4) 






GIO 


-2.7(0.4) 


81.2(5.0) 


64.5(12.4) 


-185.0(45.9) 




2ii - 2i2 


DIO 


-1.1(0.3) 


36.0 


209.6 


-94.0(28.0) 






GIO 


-0.7(0.1) 


-4.7(5.7) 


74.5(15.2) 


-58.0(14.6) 






GIO 


-0.9(0.2) 


77.9(3.2) 


34.4(7.7) 


-32.4(9.2) 


IC5179 


lio — 111 


DIO 


(1.9) 
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7 



Galaxy 



Transition 



Fif^ 



(mK) 



(km s"i) 



FWCHM/ZI)" 
(]<m s-i) 



(mK km s"!) 



NGC 7331 


lio 


- Ill 


DIO 


(0.7) 




2ii 


- 2i2 


DIO 


(0.5) 


NGC 7479 


lio 


- Ill 


DIO 


(1.1) 


IR 23365+3604 


lio 


- Ill 


DIO 


(1.0) 




2ii 


- 2i2 


DIO 


(0.3) 


Mrk331 


lio 


- Ill 


DIO 


(1.2) 



Tabic entries in parentheses are standard deviations, while entries that only list the RMS noise are non-HetRetions. 
Transitions in italies reanalyzed and/or amended with new measurements from those presented in |Mangum et al.| 



l|2008( 



Gnn / Dnn = Gaussian / Direct measurement results with nn km/s gaussian speetral smoothing. 
Heliocentrie optical velocity frame. 
^ Uncertainty in direct measurement Vhel and FWZI assumed to be 2 km/s, or 1/5 to 1/10 of a smoothed channel 
width, for all measurements. 

^ Full-width half maximum (FWHM) given for gaussian fits; full-width zero intensity (FWZI) given for direct measure- 
ments. 

® Derived from direct integration of line profile. 



ular, however, we should note that radiative excitation 
effects, in which coupling to the nonthermal continuum 
can suppress the opacity in the lowest transitions, may 
strongly affect the shape of a molecular spectral energy 
distribution (SED). To explore this effect, detailed mea- 
surements of rotatio nally-excited OH lines (see also Im- 
pellizzeri et al.|2006 ) may be worthwhile. Before obtain- 
ing interterometric nigh resolution measurements, how- 
ever, it first has to be demonstrated that the lines are 
detectable. This is shown in Figure|6] The observed OH 
line widths for the nine galaxies detected in our sam- 
ple are similar, though slightly larger, than those ob- 
tained in H2CO, suggesting a similar dynamical origin 
for the OH and H2CO emitting regions. Note, though, 
that similarity in galaxy-scale (i.e. large) line width does 
not uniquely associate the physical regions from which 
these molecular emission lines originate. Viable alter- 
nate mechanisms which can produce similarly large line 
widths include shocked gas, turbulence, and mixture of 
outflowing material with gas motions. Apparent optical 
depths for each of the OH and H2CO line region mea- 
surements range from 0.04 (Arp220) to 0.22 (IC860), 
suggesting optically thin absorption. 

Among the eight sources with detected rotationally- 
excited OH absorption, a total of five (lR01418-f 165f 
(alias niZw35), IC860, IR 15107+0724, Arp220, and 
IR 17208-0014) are weh known (Uhra) Luminous In- 
frared Galaxies ((U)LIRGs ) with prom i nent OH maser or 
megamaser emission (e.g. 



p90, ,Baanet al..l992|). The comi 



Baan 



maser emission and O 
cited ^Hi /2 level 



19891 IHenkel fc Wilson 
jmation ot ground state 



absorption in the rotationally ex- 



frared radiation field (e.g. Henkel et al. 1986 


1987 


Lock- 


ett & Elitzur 2008 


Willett et al.|201l|). We are not aware 



any 



Baan et al. 1992 1 



NGC 7331, and IR 23365+3604 (e.g. 
so that these cases are less constrained" 



4.3. Continuum Emission 

Table E lists our measured continuum levels derived 
from the mean zero-level offset of our spectroscopic mea- 
surements. Note that fo r galaxies whos e flux measure- 
ments were presented in Mangum et al. (20081 we have 
reanalyzed those flux measurements and have taken a 
mean of the baseline rather than just the central chan- 



nel offsetif] for the flux measurement. Differences be- 
tween our derived con tinuum fluxes and those quoted in 
Mangum et al. (2008) are with one exception less than 
2% for our 4.8 GHz measurements, the one exception 
being NGC 3079 where the difference is 15%. For our 
14.5 GHz continuum measurements the differences are 
less than 15%. 

For most of the galaxies in our sample t he measured 
contin u um fluxes ag r ee wit h those quoted in Araya et"aL 
( |2004| , [Baan et al."] ( |1990[ ), and the NED archive. Ex- 
ceptions to this agreement are shown in italic in Ta- 
ble [5] Continuum emission from starburst galaxies at 
these frequencies emits as a power law with spectral in- 
dex a (Sjy oc v^"). Likely emission mechanisms at the 
frequencies of our GBT observations include synchrotron 
emission (a ~ 0.7 — 0.8) and free-free emission {a ~ —2.0 
to +0.1). The GBT fluxes hsted in Table [5] suggest that 
the centimeter-wavelength emission from the majority 
of the starburst galaxies in our sample is produced via 
synchrotron emission, though several galaxies, including 
NGC 1055 and UGC 05101, appear to have flat spectra 
suggestive of free-free emission. 

5. ANALYSIS 

5.1. Comparison to Previous Measurements 

In the following we list previous measurements of the 
H2CO emission in galaxies where we have detected the 
H2CO lio — 111 or 2ii — 2i2 transitions. We also present 
a re-evaluation of our report ed detection of H2CO lio — 
111 emission in UGC 05101 ( jMangum et al.||2008[ ) which 
we now believe is a non-detection. Where available, we 
also list measurements of high spatial density molecular 

tracers in these galaxies. 

NGC 253: As noted in [Mangum et al.|(|2008], the 4.8 



and 14.5 GHz H2CO transitions have previously been de 



tected in a bsorption (Gardner & Whiteoak 1974 Baan| 
et al. 1990). Our GBT detection of the lio — lii tran- 
sition (Figure [2]) is signific antly different in shap e from 
the line profile rep orted by Gardner & Whiteoak ( 1974 ). 
Baan et al. (19971 have used the Very Large Array to 
image the 



Lio 



T^i transition, measuring an absorption 



* In [Mangum et al.||2008l l the "central channel offset" was de- 
rived by tittmg a polyriomial to the spectral baseline and capturing 
the central channel value from that baseline (continuum) fit. 
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Fig. 1. — H2CO lio — 111 spectra for galaxies toward which the 2ii 
the full-width zero intensity (FWZI) CO linewidth. 



Structure with FWHM size 39" x 12". Th e peak inte- 
grated flux measured by Baan et al. (1997) of —1.72 Jy 
beam"^ km s~^ is about 60% of the integrated flux we 
measure, suggesting that some extended structure de- 
tected in our GBT observations is missing in these in- 
terferometric measurements. P revious CS, H2CO , NH3 
and high- J CO measure ments (|Baan et a L [1990 1997| 
'Hiittcmcistcr ct al. 19 971 |Utt et al | l 2U05MGiisteii et aL 
2006||Martfn et al.|2006ffBayet et al.|2009| ) estimate spa- 



Velocity (km s ) 

-2i2 was not detected. The dotted line above each spectrum indicates 



NGC 520: |Araya et al.| ( |2004| ) first detected the H2CO 
lio ~ 111 transition m absorption, in good agreement 
with the GBT transition properties (Figure IT]). Our 
non-detection of Hllla (RMS ^ 0.25 mJy m Av 
40kms~^) is consis tent with the ^ 0.3 mJy detection of 
HllOof reported by Araya et al. (2004). Interferometric 
measurements of the CO J=l — U emission from NGC 520 



tial densities in the range lO'' cm which compares 
favorably to the io5 03±o o3 cm~^ we measure (Tabic [6|. 



by [Sanders et al.| ( |1988[ ) and |Yun fc Hibbard] ( |2001| ) re- 

veal kpc-scale structures over a region ^ 7 — —12 in 
extent. 

NGC 660: Our GBT detection of the H2CO Im - In 
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Fig. 2.— H2CO 2ii - 2i2 (panel top) and lio - In (panel bottom) spectra of NGC 253 (top left), NGC 660 (top right), Maffei 2 (bottom 
left) , and IC 342 (bottom right) . The dotted line within each spectrum indicates the FWZI CO linewidth. 
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Fig. 3. — H2CO 2ii — 2i2 (panel top) and lio — In (panel bottom) spectra of NGC 2146 (top left), the kinematical center of M82 (top 
right), M82SW, which is the (-12", -4") offset position from the kinematical center of M82 (bottom left), and NGC 3079 (bottom right). 
The dotted line within each spectrum indicates the FWZI CO linewidth. 
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Fig. 4.— H2CO 2ii - 2i2 (panel top) and lio - lii (panel bottom) spectra of NGC3628 (top left), IC860 (top right), M83 (bottom 
left), and IR 15107+0724 (bottom right). The dotted line within each spectrum indicates the FWZI CO linewidth. 
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Fig. 5. — H2CO 2ii - 2i2 (panel top) and lio - In (panel bottom) spectra of Arp220 (left) and NGC6946 (right). The dotted line 
within each spectrum indicates the FWZI CO linewidth. 



TABLE 3 
OH AND Hllla Detected Galaxies 



Galaxy 


Transition'^ 






FWZI 






(inK) 


(km s"^) 


(km s-^) 


NGC253 


Hllla 


17.6(1.7) 


219.3(3.8) 


183.5(9.0) 


IR01418+1651 


OH4750 


-5.4(0.7) 


8236.8(5.3) 


198.6(12.5) 




OH4765 


-3.2(0.6) 


8284.1(10.0) 


312.5(23.9) 


M82 


Hllla 


21.6(2.0) 


119.0(6.5) 


107.5(13.2) 






16.4(1.2) 


271.3(11.8) 


161.8(27.1) 


M82SW 


Hllla 


22.3(1.2) 


120.5(3.2) 


109.0(8.1) 






9.0(1.3) 


291.4(7.9) 


108.2(20.8) 


IC860 


OH4750 


-2.7(0.8) 


3856.0(7.1) 


132.0(17.1) 




OH 476 5 


(0.8) 






M83 


HlllQ 


2.2(0.6) 


491.8(12.3) 


135.9(30.1) 


IR 15107+0724 


OH475O 


-4.3(0.6) 


3891.4(9.6) 


190.3(24.4) 




OH 476 5 


-2.1(0.6) 


3910.5(13.2) 


203.6(35.9) 


Arp 220 


OH475O'' 


-26.3(0.7) 


5432.5(2.2) 


277.6(5.4) 




OH4765 


-11.9(0.4) 


5458.5(4.5) 


314.9(11.2) 


IR 17208-0014 


OH4750 


-3.6(0.7) 


12835.6(12.4) 


250.1(32.1) 




OH4765 


(0.7) 






IR 17468+1320 


OH4750 


-4.7(1.3) 


4737.9(7.0) 


341.5(16.7) 




OH4765 


(1.34) 






NGC7331 


OH4750 


-2.0(0.6) 


881.1(6.1) 


153.8(14.6) 




OH4765 


(0.64) 






IR 23365+3604 


OH4750 


-2.3(0.7) 


19284.6(12.7) 


398.2(35.1) 




OH4765 


(0.66) 







^ Transitions in it alics reanalyzed from measurements presented in 
[Mangum et al.H2008t . 

RMS noise levels are for 20 km s~ ^ channels. 

*^ Arp 220 OH4750 spectrum shown in |Mangum et al.lpOOS} was mis-scaled 
by a factor of 1.97. 
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Fig. 6. — OH4750 and OH4765 MHz spectra. With the exception of IR 17468 the dott ed line within eac h spectrum indicates the FWZI 
CO Unewidth. The dotted line for IR 17468 indicates the FWZI OH 1667 MHz linewidth | |Baan et al.|1993[ l. The anomalous "peaks" in the 
IR 15107 OH spectra at 3550 and 450 kms~^ are due to a fi xed-frequency receiver resonance. Note that the OH4750 spectrurn of Arp220 
shown here differs from that shown in |Mangum et aL| l |2008| due to a scaling error of a factor 1.97 in the |Mangum et al.| | |200"8| spectrum. 



transition in absorption (Figure |2]) is similar in both in- 
tensity and velocity extent as that reported by |Baan et al.| 



jMeier 



(1986), while our detection of the Ha CO 2 



Aalto et al 



IL 



^12 



tran- 



sition is the first reported 
the i^co J=l - 0, 2 - 1; 
J=l - 0; and HCN J=l - emission, while Israel (120091 



(|1995|) measured 



J=l - 0, 2-1; 



measured the ^^CO J=l - 0, 2 - 1, 3 - 2, 4 - 3 and^^CO 
J=l — 0, 2—1, and 3 — 2 emission from this galaxy. From 
CO J=l - (Aalto et al.„1995j and 3 - 2 ( |Israel|[2009 | 
images these studies derived a source diameter of ^ 15' . 



NGC891: Our 1 



Ill 



( 


2008) (Figure 




Garcia- 



1 — emission toward this edge-on galaxy. The emis- 
sion is extended along the galactic disk to spatial scales 
> 120" X 20", but is mainly concentrated in a nuclear 
condensation ^ 20" x 10" in extent. 
MafFei2: Images of the CO J=l - 0, 2 - 1, and CS 
J=2 - 1 emission (iKuno et al.||2008|), ^^CO and C^^^O 



1 - an d 2-1 and HCN J=l - emission 
et ar||2008 |), and NH3 (1,1) and (2,2) emission ([Lebron 
a 5" X 30" two- to four-component 
structure of the central bar in this nearby spiral galaxy. 
These studies indicate dense gas structures with spa- 
tial densities and kinetic temperatures ~ 10"*^ cm~'^ and 
Tk 30 — 40 K, respectively. High-excita tion NH3 
(6,6) measurements ( Mauersberger et al.||2003" l have un- 
covered a high-temperature [L'k 130 K) component 
in this galaxy. Our detections of H2CO lio — lii and 
2ii— 2i2 absorption (Figure[2| are the first reported mea- 
surements of H2CO in this galaxy. The K-doublet H2CO 
spectra are clearly composed of two velocity components 
at —98 (component "CI") and +25kms ~^ (component 
"C2"), con sistent with previous studies (Mauersberger 



et al. 



JNGO 1144: |Gao fc Solomon ( |2004a ) reported broad yet 
weak CO and HCiN J=l — emission {Tmb ~ 50 and 
2mK, respectively) toward this Seyfert galaxy. Our de- 
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TABLE 4 
OH and/or Hllla Results 



Hllla 



Galaxy 




RMS''''' 






OH 4765 
(mK) 


OH 4750 
(mK) 


HlllQ 

(mK) 


NGC 253 


3.15 


1.95 


Emis 


IC 1623 


0.91 


0.80 


0.94 


NGC 520 


0.81 


0.78 


0.74 


NGC 598 


0.68 


0.67 


0.70 


NGC 604 


1.22 


1.13 


1.23 


NGC 660 


0.79 


0.80 


0.76 


IR01418+1651 


Abs 


Abs 


0.85 


NGC 695 


0.54 


0.80 


0.55 


Mrk 1027 


0.97 


0.95 


1.02 


NGC 891 


0.94 


0.90 


0.82 


NGC 925 


0.99 


1.02 


1.07 


NGC 1022 


0.96 


0.90 


0.99 


NGC 1055 


0.80 


0.80 


0.91 


MafFei 2 


1.24 


1.12 


1.19 


NGC 1068 


2.73 


3.72 


4.22 


UGC 02369 


0.96 


0.96 


0.89 


NGC 1144 


0.78 


0.90 


0.73 


NGC 1365 


1.53 


0.55 


1.44 


IR 03359+1523 


1.17 


1.22 


1.10 


IC 342 


1.66 


1.24 


1.90 


NGC 1614 


1.14 


1.26 


1.09 


VIIZw31 


1.03 


0.71 


1.07 


NGC 2146 


2.33 


0.97 


3.18 


NGC 2623 


0.83 


0.82 


0.91 


NGC 2903 


0.42 


0.46 


0.52 


Arp 55 


1.67 


1.63 


1.68 


UGC 05101 


1.09 


1.09 


1.15 


M 82 


3.24 


3.81 


Emis 


M 82SW 


3.63 


3.63 


Emis 


NGC 3079 


2.55 


1.00 


7.61 


IR 10173+0828 


1.61 


1.68 


1.81 


NGC 3227 


2.42 


1.86 


3.28 


NGC 3627 


0.55 


0.70 


0.71 


NGC 3628 


0.94 


0.96 


0.83 


NGC 3690 


2.11 


1.87 


0.49 


Mrk 231 


1.39 


1.42 


1.26 


IC860 


Abs 


0.83 


0.90 


NGC 5194 


0.36 


0.36 


0.39 


M 83 


0.70 


0.66 


Emis 


Mrk 273 


3.63 


2.15 


4.05 


NGC 5457 


0.54 


0.70 


0.55 


IR 15107+0724 


Abs 


Abs 


0.63 


Arp 220 


Abs 


Abs 


0.91 


NGC 6240 


0.95 


0.95 


0.89 


IR 17208-0014 


Abs 


0.66 


0.81 


IR 17468+1320 


1.34 


Abs 


1.28 


NGC 6701 


0.87 


0.84 


0.60 


NGC 6921 


0.80 


0.93 


0.88 


NGC 6946 


0.55 


0.56 


0.52 


IC5179 


1.50 


1.87 


1.86 


NGC 7331 


Abs 


0.64 


0.70 


NGC 7479 


1.17 


0.87 


1.10 


IR 23365+3604 


Abs 


0.66 


0.67 


Mrk 331 


0.68 


0.72 


0.64 


^ RMS noise levels 


are for 20 km s ^ channels. 



An "Abs" entry indicates that the line was de- 
tected in absorption, "Emis" in emission (see Ta- 
ble [sj. 



tection of broad H2CO lio — In absorption (Figure [T]) 
which decomposes into three velocity components is con- 
sistent with the HCN J=l — spectra reported by |Gao fc] 



Solomon ( 2004a I . There do not appear to be any imaging 
studies of the dens e gas in this ga laxy. 
NGC 1365: Baan et al.| ( |2008[ ) report single-position 
J==l - CO, HCiN, HiNC, HCCF, N=l - CN, J=2 - 1 
CO, J=3 - 2 CS, and N=2 - 1 CN observations with 
resolutions r anging from 13" to 57" f rom th is barred spi- 



ral galaxy. Perez-Beaupuits et al. (2007) report HCN 




600 

Velocity {km s ^) 



200 400 600 800 1000 



Velocity (km s ) 



Fig. 7. — Hllla spectra. The dotted line within each spectrum 
indicates the FWZI CO linewidth. 



transitions of CO, HC N, HNC, and N: 



a-Oand N=2-l 



Sandgvist] ( |1999[ ) mapped the CO J=3 
found extended (I 



transitions of CN. Gao & Solomon (2004a) report rela- 
tively bright (160 and 8 mK, respectively) and extended 
CO and HCN J=l — emission from this galaxy, while 

2 emission and 
20" FWHM) structure dominated 
by a circumnuclear molecular torus {Otoj-us — 13"). Our 
detection of weak H2CO lio — In absorption (Figure [TJ 
is consistent with these previous molecular tracer mea- 
surements. 

IC 342: The absorption measurem ents in H2CO Im — 



111 and 2ii — 2i2 reported in Mangum et al. (20081 
were new detections of these K-doublet transitions (1* ig- 
ure[2l). Previous millimeter- wavelength H2CO measure- 



ments ( Hiittemeister et al. 1997 
Meier et al.||2011 ) estimate spatia 



Meier fc Turner||2005 

densities in the range 
10"^"" cm~'^ . High-excitation NH3 (6, 6) emission toward 
this galaxy ( Mauersberger et al.|2003 ) points to the pres- 
ence of high kinetic temp eratures. HCN J=l — im- 
ages (Downes et al. 1992) resolve IC342 into five con- 
densations with sizes ranging from 3" to 6" over a region 
- 10" X 15" in extent. 

NGC 2146: The relatively intense CO and HCN J=l-0 



emission (800 and 3 mK, respectively) measured by Gao 



& Solomon (2004a) is consistent with our detections of 
H2CO lio — 111 and 2ii — 2i2 absorption (Figure [3]), the 
first reporte d detection of this molecule in this galaxy. 
Greve et al.|(|2006|) imaged the ^^CO J=l-0 and J=2-l 



and ^CO J=l— emission from this barred spiral galaxy, 
measuring a disk-like structure ~ 5" x 20" in extent. 
UGC 05101: This was reported as a new detection of 
the li o — 111 transition in emission in [Mangum et al.| 
(20081. Additional measurements, though, have thrown 
this assignment into question as we do n ot see consistent 
results between our 2006 (reported in Mangum et al. 
)2008| ) and our new 2013 measurements. We now believe 



^12 



transitions 



J=3 — 2 emission m addition to the J=l — and J— 2 — 1 



that neither the lio — In nor the 2ii 
have been detected in this galaxy. 

M 82: The morphology of M 82's mole cular ISM is char - 
acterized by a double-lobed structure ( Mao et aL]|2000 ). 
Seaquist, Lee, & Moriarty-Schieven ( 2006p imaged the 
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TABLE 5 
Measured Continuum Levels 



Galaxy 



4.8 GHz 



14.5 GHz 





Qgrpa.b 


Arccibo^,140ft^ 


NED 


Qgrpa.b 


140ft'' NED 




( Ivl 








\^y } \^y ) 


NGC55 


PM 




0.150-0.197 


0.003 




NGC253 


1.771 


1.20 


2.0-2.5 


0.522 


0.37 


IC 1623 


0.077 




0.096 






NGC520 


0.081 


0.079 


0.087-0.126 


0.033 




NGC598 


PM 




1.300 


0.002 




NGC660 


0.151 


0.140 


0.156-0.232 


0.288 




IR01418+1651 


0.031 


0.028 




0.018 




NGC695 


0.024 




0.041 


0.011 




Mrk 1027 


0.021 










NGC 891 


0.319 




0.194-0.342 


0.101 




NGC 925 


0.030 










NGC 1022 


0.031 






0.171 




NGC 1055 


0.068 




0.063-0.086 


0.058 




Maffei 2 


0.276 




0.243-0.375 






NGC 1068 


1.680 




1.770-2.187 


0.555 


0.680 


UGC 02369 


0.112 










NGC 1144 


0.062 




0.052 


0.011 




NGC 1365 


0.235 




0.180-0.230 


0.052 




IR 03359+1523 


0.016 










IC 342 


0.086 




0.080-0.277 


0.050 




NGC 1614 


0.072 




0.063 






vnzwsi 


0.032 






0.015 




NGC 2146 


0.389 






0.223 


0.30 


NGC 2623 


0.061 


0.064 


0.057-0.073 


0.009 




Arp55 


0.039 










NGC 2903 


0.115 


0.041 


0.118-0.180 


0.019 




UGC 05101 


0.054 




0.077 


0.048 




M82 


3.233 


3.55 


3.67-3.96 


1.304 


1.79 1.73 


M82SW 








1.128 




NGC 3079 


0.342 


0.33 


0.32 


0.159 


0.14 


IR 10173+0828 


PM 


0.020 








NGC 3227 


0.559 


0.008 


0.045 






NGC 3627 


0.121 










NGC 3628 


0.186 


0.131,0.33 


0.222-0.283 


0.043 


0.11 


NGC 3690 


0.321 




0.362-0.407 


0.115 




NGC 4631 






0.299-0.438 


0.033 




NGC 4736 






0.090-0.119 


0.022 




Mrk 231 


0.321 




0.417 


0.129 




NGC 5005 






0.060 


0.001 




IC860 


0.010 




0.042 


0.031 




NGC 5194 


0.150 




0.215-0.436 


0.023 




M83 


0.322 




0.648-0.780 


0.071 




Mrk 273 


0.012 




0.070-0.103 


0.028 




NGC 5457 


0.023 






0.013 




IR 15107+0724 


0.017 


0.031 




0.021 




Arp220 


0.217 


0.172,0.22 


0.22 


0.082 


0.11 


NGC 6240 


0.141 




0.178 


0.035 




IR 17208-0014 


0.071 


0.062 


0.068 


0.028 




IR 17468+1320 


0.190 




0.047 






NGC 6701 


0.025 




0.022 


0.008 




NGC 6921 


0.028 




0.027 


0.035 




NGC 6946 


0.167 




0.219-0.531 


0.030 




IC5179 


PM 




0.079 






NGC 7331 


0.095 


0.040 


0.070-0.096 


0.007 




NGC 7479 


0.038 




0.040 






IR 23365+3604 


0.015 




0.010 


0.001 




Mrk 331 


0.028 


0.025 


0.028 







^ GBT mea s urements shown in italic s differ from the measurement range defined by 
[Araya et al.| ) |'2004[ l, [Baan et al.| ) |l990[ l, and NED measurements by more than 25%. 

PM = Poor measurement due to insufficient total power stabihty. 

The Arecibo 4.8 GHz continuum levels were obtained from spectral baselines by |Araya| 
|et al.| | [2004l l, their Table [l] 

The NRAO 140' 4.8 and 14.5 GHz continuum levels from iBaan et al.l l|l990l, their 
Table [1] ' ' ^ ' 
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CO J=:6 — 5 emission from this galaxy, finding emis- 
sion extended over a disk-like structure 40" x 15" in ex- 
tent. H2CO was p reviously detected in a bsorption in the 
111 line by 



lio 

the 2 ii 



and 



Graham et al. (|1978|), in emiss ion in 

^03 



2o2 ('JrsTTEE) Hues byjBaan 
et al.|(|1990|), and in emission in the 2o2 — loi (146 UMzj, 



2^7and 822 - 221 / 821 - 220 (218 GHz) lines by 



Miihle et al. (l2007|). We detect both the 4.8 and 14.5 
lbs 



GHz lines in absorption with high confidence (Figure ^ . 
Investigation of our LVG model predictions of the rel- 
ative intensities of the 4.8 and 14.5 GHz emission (see 
^5.8|) cannot reproduc e the 14.5 GHz emission observed 
byTBaan et al.| (|1990|). The HllOa and HI 11a lines, as 



well as other higher frequency hydrogen radio recombi 



nation lines, have been detected in M82 (iSeaguist, Bell, 



fc BignelllflMl |Bell et al.| [19841 |Rodrigu cz-Rico et a 



20041 



We re-detect the Hill 



two components identified by Bell et al 



ine and co nfirm the 
([T9'84|). Pre- 



vious CS and H2CO mea surements ( Baan et al. 



Hiittemeister et al 
the range 10^ 



1990 



19971 estimate spatial densities m 



cm 

Recent Herschel-SPIRE imaging spectroscopic mea- 
surements of the molecular emission toward M82 



(Kamenetzky et al. 2012) have characterized the phys- 
ical conditions in this galaxy using measurements of 
^^CO and -"^^CO emission over a wide range of molec- 
ular excitation. By using a LVG radiative transfer 
model with Bayesian likelihood analysis of the J=4 — 8 
throu gh J=18 — 12 ^^CO and ^'^CO emission toward 



M 82 |Kamenetzky et al. (20121 confirmed the existence 



of two temperature components (Panuzzo et al. 2010): 
A cool component with median i'x = 85K [1-a range 
12-885 K) and n(H2) = lO^-^* cm-^ (1-cr range lO^-^S- 
2^05.15 Q]-|2-3^ g^jjjj ^ warm component with median Tk = 
486 K (1-cr range 844-548 K) and n(H2) = lO^-^^ cm"^ (l- 



a range 10 



3.17_ 



10 



3.96 



cm ■^). Our H2CO measurements 



clearly correspond to the cool component measured by 
Kamenetzky et al. (2012), with our derived spatial den- 
sityofnpl^y^TO^-''^^^® cm~^ overlapping at the high- 



density limit of their derived 1-a range. 

We also measure a position offset of (—12,-4) arcsec 
from the nominal M82 position in the H2CO 2ii — 2i2 
transition only, which we refer to as M82SW (Figure Isl). 
This is the "P8" HON J=l - peak position noted by 



Brouillet fc Schilke ( 1998 1, and corresponds to within 2" 
of the "southwestern molecular lob e" where NH3 emis - 
sion was detected in this galaxy (Weifi et al. 2001). 
M 82SW appears to be ^ 50% stronger in H2GO emission 
than the nominal galaxy center position (Table |2| , but 
note that the spatial resolution of our H2CO measure- 
ments {9b = 158" and 51") allows for some sampling of 
both components in all of our M 82 measurements. I n the 



NH3 measurements of M82 ( Mangum et ar] 2018) only 
the M 82SW position possesses detectable iNHs emission, 
and is the source of the kinetic temperature assumed in 
our LVG analysis (j |5.8 ). 

NGC 3079 : Reported as a nondetection in Mangum 
et al. (20081, additional integration resulted in a detec- 
tion of the lio — 111 Ihic (Figur e [3|, consiste nt with the 
absorption line reported by Baan eFal. ( 1986( —2.1 mJy) 



[Pirez-Beaupuits et al.) ( |2d07p detect CU J =1 - and 
J=2 - 1, CN M=l - aSa^ - 1, and HON and HNC 
J=l — and J=3 — 2 emission from the nucleus of this 



galaxy. From the HCN J^l — image pr e sented by 



Kohno et al.| ( |2000[ ) Perez-Beaupuits et al.| ( 2007|) de 
rive a sour ce size ot 5 x 5". 450 /zm, 850 /im ( [Stevens fc 
Gea^[2000] ), and 1200 /im continuum imaging reveals an 
unresolved core ~ 20" in extent. 

NGC 3079 is one of two galaxies in our sample which 
shows complex absorption and emission structure in its 
H2CO lio — 111 spectrum. The lio — In and 2ii — 2i2 
transitions are dominated by velocity components at 
^ 1010 and ~ 1115 km s~^ (which is the systemic ve- 
locity), with the lio — In transition exhibiting a tenta- 
tive third broad component at ~ 1500 km s~^. The two 
velocity components which are apparent in both H2 CO 
transitions are de tected in CH3OH, OH, and HI (Impel- 
lizzeri et al.|2008 ). The excitation characteristics ot these 
three velocity components differ dramatically: 

• lOlOkms^^ Component (CI): lio — In emission 
with 2ii — 2i2 absorption suggests high spatial den- 
sity. In subsequent discussion we refer to this as 
component 1 (CI). 

• lllSkms^^ Component (C2): Complex lio — In 
emission and absorption spectral structure suggest- 
ing that this component has a spatial density and 
kinetic temperature which lies near the point where 
the lio — In transition goes from absorption to 
emission. This lio — In spectral structure with 
2ii — 2i2 absorption suggests lower spatial density 
than the 1010 km s^^ component. In subsequent 
discussion we refer to this as component 2 (C2). 

• 1500 km s^^ Component (C3): Only detected in 
broad lio — In absorption at 4.5cr confidence in 
peak intensity. In subsequent discussion we re- 
fer to this as component 8 (C3). We are not 
aware of corresponding spectral line detections of 
the 1500 km s""'^ component in other high density 
molecular tracers (though the CO emission toward 
NGC 8079 includes this velocity; Figure |3]). 



1 



NGC 3628: The 

absorption b y [Baa n et al. 
Araya et al. ('2004) 7wEii rt' 



I n tra nsition was detected in 
1990) and confirmed by 
— 2i2 transition was not 
detected byiBaan etTaTl ([19901). The GBT lio - In line 



( fl986l 

3ne 2ii 



prop erties show good a greement with those measured by 



both 



(200. 



Baan et al 



(1986 S = -2.5 mJy) and Araya et al 



97 ± 0.21 mJy) to withi n the calibra- 
tion uncertainties (Figure [4|). Israel ( 2009[ ) measured the 
i^CO J=l - 0, J=2 - 1, J=3 - 2, J=4 - 8 and ^^CO 
J=l — 0, J=2 — 1, and J=8 — 2 emission from this galaxy. 
From their CO J=3 — 2 map one notes that the source 
size is < 15". 450 u m and 850 /im continuum imaging 
([Stevens e t al.[[2005 ) expose a '-^ 80" core. 

I ll emission l ine w as detected by 



860: THe T 



Baan et al. ( [1998[ ) and [Araya et al.| ( |2004[ ), and the GBT 

detection snows good agreement m line properties (Fig- 
ure |4|. Our discovery of H2CO 2ii — 2i2 absorption is a 
new result. 15" res olution CO J=2 — 1 and 1-0 emis- 
sion measurements |Yao & Seaquist ( 2008 ) provide an 



estimate of the maximum source size 
NGC 5194 (M51): The "Whirlpool Galaxy" is a well- 
studied Sbc galaxy with a weak AGN. The inner 2. '5 of 
this galaxy have been imaged in a variety of mol ecular 
tracers, including CO, "CO, and C^O J=l - ( ,Kohno 
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iFaD (11996]);^ 

et al. 2010 T, and HCJN an5Tir !0+ J=l - (jK'ohno et al. 



Aalto et al. (1999)), CO J=2-l (Schinnerer 



( |T996| ); |Scfiinnerer et al. ( 2010| ). These dense gas studies 



have revealed structure associated with the nearly face 
on spiral arms down to spatial scales ^ 120 pc. Our 
H2CO lio — 111 absorption measurement (Figure [T]) is 
the first reported detection of H2CO in this galaxy. 
M83: Both H2CO transitions are detected in absorp- 
tion, with the 14.5 GHz transition showing a significantly 
larger optical depth than the 4.8 GHz transition (Fig- 
ure |4]). CO J=2 — 1 an d 3 — 2 emission imagi ng of this 
barred starburst galaxy ( Sakamoto et al.|2004 | indicates 
a 45" X 15" structure. 

Mrk273: Ima ged in CO J=3-2 by|Wilson et a l.|( |2008 ) 



Based on th ese Wilson et al. ( 2008 ) measurements |lono 
et al. (20091 deriv e a deconvolved CO J 



source 



0.3". |Gracia-Carpio et al.| (|2008 ) measured 

IT and 3 -2, 
Fur- 



size of 

moderately intense MCiN and HCO"*" J 
along with CO J=l — emission from this ULIRG 
thermore, Herschel has detected OH emission due to out 



flows from this galaxy ( Fischer et al.|[2010 ). Our H2CO 
lio ~ 111 absorption measurement is the first reported 
measurement of H2CO in this galaxy (Figure [l|. 
IR 15107+0724: Our Im — In measurement is consis 



tent w ith th at reported by Baan et aL] ( 1993 1 and Araya 

{2004| . We also detect the 2 ii — 2 i2 line in absorp 



gureW|. Planesas et al. (1991) imaged the CO 



et al. 

tion ( 

J=:l — emission from this galaxy, deriving a compact 
nuclear source with size ~ 3". 
Arp220/IC4553: The H2CO Im 



detected in emission by Baan et al. ( 1986 



transition was 
but the 



2ii — 2i2 line was not subsequently detected ( |l^aan et al. 



1990). We find good agreement between our GB'l' mea- 



surement (Figure [5 ) and the total integrated intensity 
of the lio ~ 111 sub-arcsecond resolution interferomet- 



ric measurements reported by Baan & Haschick ( 1995 ) 



Our 4.8 and 14.5 GHz GBT spectra ar e consistent with 



those reported by Araya et al. J 2004 ). Note also that 
previous studies (Baan et al. T^Sb^ have suggested that 
H2CO lio — 111 emissi on is produced b y maser amplifi- 
cation in this galaxy. In Mangum et al. (2008) we argued 



that, to the contrary, the lio ~ lii emission in Arp220 
is simply the signature of a dense gas component in this 
ULIRG. 



Mangum et al. (2008) summarized the previous molec- 
ular spectral line measurements of Arp 220. Of partic- 
ular note is the dichotomy between the relative inten- 
sities of various dense gas tracers from the two nuclei 
which comp rise this merging system. As noted by |Greve| 



et al. (2009), this dichotomy reflects differing spatial den- 
sities m these two components: the western component 
possessing lower spatial densities, the eastern component 
possessing higher spatial densities. Our H2CO measure- 
ments partially reflect this trend. Both H2CO transi- 
tions can be flt with velocity components at ^ 5330 
and ~ 5460 km s^^, likely corresponding to the west- 
ern and eastern velocity components noted previously. 
We also find an indication of a third component at 
^ 5600 km s~^. We are not aware of corresponding 
spectral line detections of this ^ 5600 km s~^ velocity 
component in other molecular species. 
Recent Herschel-SPIRE imaging spectroscopic m easure- 



ments of the mo lecular emission toward Arp 220 ( Rang- 
wala et al.|| 2011| ) have characterized the physical con- 



ditions in this galaxy using measurements of ^^CO and 
HCN emission over a wide range of molecular excita- 
tion. By using a radiative transfer model with Bayesian 
likelihood analysis of the J=4 — 3 through J=1 3 — 12 



^^CO emission toward Arp 220 Rangwala et al. (2011) 
derived two temperature components! A cool compo- 
nent with median Tk = 50K (1-ct range 34-67 K) and 



2.6 



n(H2) = 10-^**cm"'^ (1-cr range 10 
a warm component with median Tx 
range 1247-1624 K) and n(H2) = lO^ ^ 



10^-^ cm ^) and 
= 1343 K (1-cr 
cm"*^ (1-(T range 
10'^ °-10'^-^ cm~'^). A similar LVG radiative transfer 
and Bayesian likelihood analysis of the HCN J=12 — 11 
through J=17— 16 emission yielded a single-temperature 
fit with Tk > 320K and n(H2) > lO^-^cm-^. Our 
H2C0-derived spatial density (n(H2) = W^-O^^om ^^-3-^ 
and NHs-derived kinetic temperature (T^ = 234 ± 52 K) 
correspond most closely to the s patial density and kinetic 



temperature derived from the Rangwala et al. ( 2011| ) 



HCN measurements, though are slightly lower in both 
quantities (see Table [6t . The cool and w arm components 
derived from the |Rangwala et al. (2011) CO analysis do 
not correspond to any physical components in our H2CO 
or NH3 measurements, suggesting that the contribution 
to the CO emitting gas in Arp 220 from high spatial den- 
sity components is small. 

NGC6240: Our GBT detection of the lio - In tran- 
sition in absorption (Figure [ij 5.2a in a single smoothed 
20 km s~^ c hannel) is not cons istent with the emission 



et al. 



19931 



As noted by Mangum 
11 spectrum peaks near 730U 



reported by Baan et al 
' (f2008^, our lio 

^ slightly blues hifted relative to th e systemic ve 



our 

km s 

locity o f 7359 km s"^ (|Gr eye et al [ | 2009[) . CO J=l - 
images ( Bryant fc Scoville|l999| ) indicate a molecular gas 
source size ot lU". NGC 6240 has also been imaged in 
CO J=3-2 and HCO+ J=4 - 3 by [Wilson et al.| (2008 ). 

(I2OO8I) 



Based on the Wilson et al. 



lono 



measurements 

et al. (2009) derive a deconvolved CO J=3 — 2 source 
slieof 0:^ 0.6". The HCO+ J=4 - 3 image presented 



by Wilson et al. (2008) suggests a source size < 0.5" 



NGC 6921: Gao & Solomon 
HCN 



( 2004a ) measure very weak 
and 1 mK, respectively) 



CO and HCN J=l - (- 
emission toward this spiral galaxy. We are not aware of 
any spectral line imaging measurements of this galaxy at 
radio wavelengths. Somewhat surprisingly, given the low 
intensity of the HCN J=l — measurement, we detect 
significant H2CO lio ~ In absorption toward this galaxy 
(FigurelT]). As the 4.8 GHz continuum flux is measured to 
be only z8 mJy toward this galaxy, H2CO absorption of a 
strong background continuum source does not provide a 
mechanism for producing significant H2CO absorption. 
The H2CO lio — 111 absorption we measure is clearly 
absorbing the co smic microwave backg round. 
NGC 6946: In [Mangum et al.| ( |2008| we reported the 
discovery of the H2CO lio — In transition in absorp- 
tion. We add to this the H2 CO 2ii — 2i2 transition, 
also in absorption (Figure Is l). 
^ (pOOTF i 



and 
J=2 



Schinnerer et al 
1, and HCN J 



Schinnerer et al . 
imaged the CO J 



(^2006) 



1 — (J emission from this galaxy, 
detecting a compact nuclear source size of ~ 2" and a 
"nuclear spiral" structure 5" x 10" in size. 
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The apparent peak op tical depths (r, see Equation 2 
of Mangum et aL 20081 for our H2CO absorption mea- 
surements ('l'abie[2|) have been calculated using the GBT 
continuum emission intensities listed in Table [5] and as- 
suming Tex < Tc. With the exception of NGC6921, 
which has an optical depth of ~ 0.06, all of the other 
galaxies detected in H2CO in our sample have optical 
depths < 0.007, indicating optically-thin H2CO emission. 

5.3. Spatial Density and Column Density Derivation 
Using LVG Models 



As was done in Mangum et al. (2008), to derive the 
H2 spatial density (number density) and H2CO column 
density of the dense gas in our galaxy sample, we us e a 
Large Velocity Gradient (LVG) model ( |Sobolev|1960l ) to 
calculate the radiative transfer properties ot the H2CO 
transitions. The detailed properties of our im plementa- 
tion of the LVG ap proximation are described in |Mangum| 
& Wootten ( 1993 ) . One important point regarding this 
implementation of the LVG model is our scaling of the 
calculated ortho-H2CO/He excitation rates to those ap- 
propriate for collisions wit h H2. Following the recom- 
mendation of Green (1991), we scale the calculated He 
rates by a factor ot 2.2 to account for (1) the reduced 
collision velocity of He relative to H2 , which scales as the 
inverse-square-root of the masses of He and H2 an d (2) 
the la rger cross section of the H2 molecule 1.6; Nerf 
(|1975|)) relative to He 



Radiative transfer models of the molecular emission 
in astrophysical environments are dependent upon colli- 
sional excita tion rates for the mol e cule(s ) un der study. 



siona. 



As no t ed by Mangum & Wootten (1993) and Mangum 
et al.| (|2008|) , the uncertainty associated with the colli- 



excitation rates must be considered in any analysis 
of the physica l conditions d erived from radiative trans- 
fer modelling. Green ( 1991 ) suggests that the total col- 
lisional excitation rate tor a given H2CO transition is 
accurate to ~ 20%. This implies that the physical condi- 
tions derived from our LVG modelling are limited to an 
accuracy of no better than 20%. 

For the 13 galaxies and one galaxy offset position 
(M82SW) where both the H2GO lio - In and 2n - 2i2 
transitions were detected we can derive a unique solu- 
tion to the (ensemble average) spatial density and H2CO 
column density for an assumed gas kinetic temperature. 
This unique solution to the physical conditions is derived 
by fitting to the intercept between the H2CO lio — In 
and 2ii — 2i2 transition ratio and the H2CO lio — In 
transition intensity at the assumed kinetic temperature. 
Absorption line measurements afford the possibility of 
using the measured apparent optical depths in the LVG 
model fitting procedure. For our galaxy sample, where 
with the exception of M 82 Tc < Tcmb, this method has 
the disadvantage of involving an estimate of the excita- 
tion temperature (Te^) in order to derive the apparent 
optical depth. Since T^x is dependent upon the input 
physical conditions that we are attempting to derive, we 
do not believe that fitting to measured apparent H2CO 
optical d epths is advantageous . In fact, as we showed in 
§ 5.3.3 of Mangum et al. (2008), these two approaches to 
estimating the excitation ot H2CO absorption lines in the 
presence of the cosmic background radiation and weak 
ambient continuum emission are roughly equivalent for 
the measurements presented in this paper. Furthermore, 



we showed in Mangum et al. (2008) that for M82, which 
possesses strong continuum emission, ignoring contribu- 
tions due to background continuum emission produces a 
derived spatial density which is a lower-limit to the true 
average spatial densi ty in this galaxy. 



As was done injMangum et al.| (|2008|), a model 



grid of predicted H2GU transition intensities over a 
range and step size in spatial density, ortho-H2CO col- 
umn density per velocity gradient, and kinetic tem- 

-l4.0 



perature (n(H2), N(ortho-H2CO)/Ai;, Tk) 



^10.0 



10 



14.0 



cm 



-2 



/km s 



-1 



20- 



= (10^ 
300 K) and 



lO'-^cm--", 10 

(Alog(n(H2)), Alog(N(ortho-H2CO)/Aw), ATk) = 
(0.03 (cm-3), 0.04(cm-Vkms-i), 5K) was calculated. 
Note that the model grid used in the present calcula- 
tions is significantly la rger in all three mod eled parame- 
ters than that used in Mangum et al. ( |2008| ). The pre- 
dicted transition intensities were then compared to our 
measured HtCO Iiq — In and 2n — 2i2 transition inten- 
sities (Table [2]). 

Table [6] lists the derived LVG model best-fit physical 
conditions for all of the H2C0-detected galaxies in our 
sample. There are three different sets of assumptions 
used to derive these physical conditions dependent upon 
the information available: 

• Both Transitions Detected: For the 13 galaxies and 
one galaxy offset position where we have detected 
both the H2CO lio ~ In and 2ii — 2i2 transitions, 
we calculate explicit values for the spatial density 
and ortho-H2CO column density assuming Tk — 
Tdust, with Tdust from Table[T] For the 12 galaxies 
where we also have NHa-baseq kinet ic temperature 
measurements (Mangum et al.|2013D , we also derive 
spatial density and column density assuming these 
gas-derived kinetic temperatures. As NGC 253 and 
IC 342 possess two NHs-derived temperature com- 
ponents LVG model fits assuming both of these ki- 
netic temperatures are listed. For the LVG model 
fits with NHs-derived kinetic temp eratures, kinetic 



temperature uncertainties (Mangu m et al.||2013 l 
are included. By using NHs-derived kinetic tem- 
peratures to constrain our H2C0-derived spatial 
densities we are assuming cospatiality for these two 
dense molecular gas tracers. As the spectral line 
profiles from NH3 and H2C O are similar in all but 
one galaxy (NGC 660; see S 5.3.1), the assumption 
of cospatiality seems plausible. 

• One Transition Detected and Limit to Other Tran- 
sition: For the 8 galaxies and one galaxy velocity 
component where we have detected only one of the 
two K-doublet transitions, but also have a limit to 
the nondetected transition, we make two estimates 
of the physical conditions: 

— We derive limits to both the spatial den- 
sity and ortho-H2CO column density assum- 
ing a kinetic temperature of 40 K, a 3ct in- 
tensity limit, and an equivalent line width 
for the nondetected transition. Furthermore, 
since our LVG modeling results imply that the 
lio ~ 111 transition goes into emission before 
the 2ii — 2i2 transition as the spatial density 
is increased (for a fixed kinetic temperature). 
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a detection of absorption in the lio — lii tran- 
sition means that the (undetected) 2ii — 2i2 
transition must also be in absorption. On the 
other hand, a detection of the lio — In transi- 
tion in emission implies that the 2ii — 2i2 tran- 
sition can either be in emission or absorption. 
These results are listed as limits in columns 3 
and 4 of Table 1 



— As was done in |Mangum et al.| ( |2008[ ), we cal- 
culate the LVG model derived column density 
assuming n(H2) = 10^ cm~^ (for absorption 
sources) or 10^-^^ cm~'^ (for emission sources). 
These assumed values for the spatial den- 
sity are reasonable averages from our H2CO 
galaxy sample where spatial density was de- 
rived. We also assum^ Tk — Tdust ■ This 
line of reasoning allows us to derive an ac- 
tual value rather than a limit for N(ortho- 
H2C0)/Av that we can use in a calculation 
of the dense gas mass. The column densities 
derived using these assumptions are listed in 
column 3 of Table |6l 

• Neither Transition Detected: For the 35 galaxies 
where we derive only limits to the lio — In and/or 
the 2ii — 2i2 transition intensities we assume a 
3(7 limit to the measured yet undetected intensi- 
ties. This allows us to calculate a limit to the 
ortho-H2CO column density assuming n(H2) = 10^ 
cm~'^ and Tk ~ T^dust- Note that we conservatively 
do not assume a line width in these limits, which 
would decrease the intensity limit by the square- 
root of the number of channels over which the line 
is integrated. 

For all of our LVG calculations we assume Tcmb — 2.73 K 
and negligible contribution due to any backgroun d con- 
tinuum emission (T^, = 0). In § 5.3.3 of ,Mangum et al. 
([2008) we showed that our results represent a lower-limit 
to the true spatial density when contributions due to 
strong background continuum sources (S > 1 Jy) are ne- 
glected. 

5.3.1. The Multiple Dense Gas Components of NGC 660 

Comparison of the NH3 and H2CO K-doublet spectra 
toward NGC 660 suggests that the bulk of the emission 
from these two molecules originates from different spatial 
components. The NH3 (J,K) = (1,1) through (7,7) (for 
J=K) and (2,1) rotation-inversion transitions have been 
detected toward NGC 660 with multiple absorption com- 
ponents comprising a total FWZI ^ 50—100 km s~ ^ cen 



tered near V/ie; 840 kms^^ (Mang um et al.|2013| ). The 
Vhei and FWZI values determined from the IN H3 spec- 
tra correspond to the narrow absorption peak embedded 
within the H2CO 2ii — 2i2 absorption component (Fig- 
ure l2| . Analysis of the (narro w) NH3 emission within 
NGC660 (Mangu m et al.j2013[ ) strongly suggests that it 
contains four velocity components, two of which originate 
in the disk of the galaxy. The bulk of the H2CO absorp- 
tion traces the nuclear region of this galaxy. In spite 



of the partial non-cospatiality of the NII3 and II2CO ab- 
sorption in this galaxy, we present physical condition cal- 
culations assuming both the dust- and NH3-derived ki- 
netic temperature for H2CO components associated with 
the disk component. 

5.3.2. LVG Model Dependence on Kinetic Temperature 

It is difficult to determine kinetic temperatures in ex- 
tragalactic molecular clouds. The most common kinetic 
temperature probe is interstellar dust, which is available 
for a large number of galaxies. However, dust and gas 
kinetic temperatures, even though being about 10 K in 
local dark clouds, are not always the same. Good cou- 
pling is normally only achieved at densities in excess of 
10^ cm^'^. In the case of substanti al gas heating through 
cosmic rays ( Papadopoulos 20101, even high densities 
cannot prevent a signihcant discrepancy between dust 
and gas temperatures. 

Ammonia (NH3) is a more direct probe of the gas. Its 
temperature sensitive inversion transitions trace molec- 
ular gas in star forming as well as in quiescent region s 



(e.g., Benson & Myers 1983 Mauersberger et al. 2003) 



As is the case with many molecules, INH3 shows some 
chemical peculiarities in star formation regions in our 
own Galaxy, which have to be kept in mind. NII3 is easily 
destroyed in photon dominated regions and shows partic- 
ularly high abundances in hot cores, where temperatures 
of 100 K or more lead to dust gr ain mantle evaporation 



(e.g., Mauersberger et al. 1987). Differences in abun 



dances in individual Galactic regions can amount to two 
orders of magnitude. This is also the case when compar- 
ing NII3 abundances between M82 and other detected 
galaxies. Unlike dust, ammonia has been measured in 
only a limited number of galaxies. 

Formal dehyde is another direct pr obe of kinetic tem- 
perature (Mangum & Wootten||1993). T he first applica- 
tion o f the technique described oy M angum fc Wootten 
( |1993D , using the H2CO 3o3-2o2, 322-221, and 321 —220 
transitions, has been applied to the kinetic temperature 
measurement of M 82 by Miihlc et al. (20071). 

As noted in Mangum et al.| (|2008), our LVG model- 
based derivations of the spatial density and ortho-Il2CO 
column density are dependent upon the assumed kinetic 
temperature used in these models. As we did not have 
a complete set of gas kinetic temperature measurements 
for th e samp le of starburst galaxies presented in |Mangum| 
et al. (2008), we assumed Tk = Tdust, which were gen- 
erally in the range 30- 40 K. An update of Figure 10 from 
Mangum et al. (2008), showing a graphical representa- 
tion of our LVG model derivations of the spatial den- 
sity and ortho-H2CO column density assumi ng Tj^ = 
Tdust is sho wn in Figure |8] These results of [Mangum] 
et al.| ('2008) naturally leaa to the conclusion that spa- 
tial density is an important factor in the development 
and evolution of starbursts in galaxies. The higher mean 
spatial densities of > 10^-^ cm~^ found in the LIRG 
and ULIRG (Arp22q , IR 15107-^0724) measurements of 
Mangum et al. ( 2008 ) contrast with the lower mean spa- 
tial densities of lO'* — 10^'^ cm~'^ found in the normal 
starburst galaxies of . 



As we now have IN 



(2008) 



5 We assumed Tk = 40K in [Mangum et al.| | [2008| . Note, 
though, that T^^^t ~ 40 K for most galaxies in our sample. 



' Mangum et al. 

JM3 measurements of the dense gas 
kinetic temperature in 11 galaxies with a total of 14 ve- 
locity components (iMangum et al. 2013[ ), sampled over 
similar though slightly smaller [Ub — 30 ) spatial scales 
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TABLE 6 






Derived Physical Conditions 








log(N(ortho-H2CO)/Av)''''^ 


iog{n[^rL2 ) ) 




(K) 


(cm~^/km s~^) 


(cm ) 


NGC 55 


Td„,t = 40'' 


< 11.31 




NGC 253 


Td„et = 34 


12.14 ± 0.04 


5.11 ± 0.02 




Tjf (NH3) = 78 ± 22 


12.24 ± 0.06 


5.03 ± 0.06 




TjcCNHa) > 150 


> 12.54 


< 4.82 


IC 1623 


Td„et = 39 


< 10.71 




NGC 520 


Tdust = 38 


< 10.84, 11.13 ± 0.08 


< 4.58 


NGC 598 


Tdust = 40'' 


< 10.87 




NGC 660 


Td„et = 37 


11.53 ± 0.16 


5.34 ± 0.04 




Tj<(NH3) = 160 ± 96" 


11.94 ± 0.48 


4.88 ± 0.48 


IR01418+1651 


Td„,t = 40'' 


< 10.72 




NGC 695 


Td„et = 34 


< 10.88 




Mrk 1027 


Tdust = 37 


< 11.07 




NGC 891 


Tdust = 28 


< 10.66, 10.90 ± U.lu 


< 4.50 


NGC 925 


Tdust = 40'' 


< 11.12 




NGC 1022 


Tdust = 39 


< 11.05 




NGC 1055 


Tdust = 29 


< 10.88 




Maffci 2 


Tdust = 40'' 


11.17 ± 0.12 


4.93 ± 0.09 


Maffci 2C1 


Tk(NH3) = 62 ± 25 


11.27 ±0.87 


5.00 ± 0.50 


Maffci 2C2 


Tjf (NH3) = 64 ± 24 


11.00 ± 0.50 


4.91 ± 0.27 


NGC 1068 


Tdust = 40 


< 11.00 




UGC 02369 


Tdust = 40'' 


< 11.00 




NGC 1144 


Tdust = 40'' 


< 10.91, 11.31 ± 0.02 


< 4.16 


NGC 1365 


Tdust = 32 


< 10.71, 10.83 ± 0.06 


< 4.85 




Tdust = 40" 


< 11.08 




IC 342 


Tdust = 30 


11.27 ± 0.22 


5.05 ± 0.11 




Tjf(NH3) = 24± 7 


11.33 ± 0.28 


5.02 ± 0.14 




Tjf(NH3) > 150 


> 11.37 


< 4.88 


NGC 1614 


Tdust = 46 


< 11.16 




VIIZw31 


Tdust = 34 


< 11.06 




NGC 2146 


Tdust = 38 


11.25 ± 0.12 


5.33 ± 0.03 


NGC 2623 


Tdust = 40" 


< 10.97 




Arp 55 


Tdust = 36 


< 12.42 




NGC 9Q03 


Tdust = 29 


< 10.59 




TTGC 5101 


Tdust = 36 


< 10.89 




M 82 


Tdust = 45 


11.80 ± 0.02 


4 95 + 02 




Tk(NH3) = 58 ± 19 


11.84 ± 0.06 


4 09 _|_ n nc 


M 82SW 


Tdust = 45 


11.86 ± 0.08 


c no _|_ f) no 




Tjf (NH3) = 58 ± 19 


11.86 ± 0.08 


5.05 ± 0.08 


NGC 3079C1 


Tdust = 32 


12.34 ± 0.28 


5.56 ± 0.02 




Ta'(NH3) > 100 


> 12.46 


< 5.39 


NGC 3079C2 


Tdust = 32 


11.83 ± 0.10 


5.46 ± 0.02 




Tif(NH3) > 150 


> 12.22 


< 5.09 


NGC 3079C3 


Tdust = 32 


< 11.07, 10.48 ± 0.08 


< 5.33 


IR 10173+0828 


Tdust = 40" 


< 12.36 




NGC 3227 


Tdust = 40" 


< 12.33 




NGC 3627 


Tdust = 30 


< 10.90 




N GO 3628 


Tdust = 30 


11.19 ± 0.10 


4. DO ± U.i4 


NGC 3690 


Tdust = 40" 


< 10.86 




NGC 4631 


Tdust = 30 


< 10.92 




NGC 4736 


Tdust = 40" 


< 11.33 




Mrk 231 


Tdust = 47 


< 11.02 





to our H2CO measurements, we have used these values 
to constrain our LVG models and derive revised spa- 
tial densities and ortho-H2CO column densities. Fig- 
ure [9] shows the LVG model-predicted spatial densities 
and ortho-H2CO column densities for these 11 galaxies 
(and their individual velocity components, as appropri- 
ate). In all 11 galaxies with a total of 14 velocity compo- 
nents the measured kinetic temperature is significantly 
higher than the previously-assumed Tdust, thus driving 
the derived spatial densities significantly lower. The cor- 
relation between higher kinetic temperature and lower 
spatial density is shown in Figure [TOj The range of best- 
fit values for spatial density and column density for each 
galaxy listed in Table [6] and shown in Figure [9] are driven 
by signal-to-noise limitations of our H2CO measurements 
and the uncertainties in our derived kinetic temperature 
measurements . 



Mangum et al. C2008') note, based on LVG model- 
derived spatial densities which assumed Tk — Tdust, that 
there appeared to be a correlation between infrared lumi- 
nosity and spatial density. This was purported to b e an- 
other representatio n of the Lm—Mdense correlation ( Gao 
& Solomon'2004b). The revised spatial density measure 



ments, which include appropriate dense gas kinetic tem- 
perature measurements, now point to a relatively narrow 
range in spatial density of lO"*^'^ — 10^'^ cm~^ in our star- 
burst galaxy sample. We should stress that while the 
assumed kinetic temperature influences the spatial den- 
sity that we derive from our H2CO measurements, the 
density-kinetic temperature anti-correlation is dramati- 
cally smaller in parameter space than most other molec- 
ular tracers; only 1.2 dex in spatial density (Figure |9|. 
In ^|6]we analyze the implications of the narrow range m 
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TABLE 6 — Continued 



Galaxy 




log(N(ortho-H2CO) /Av)^'" 


log(n(H2))'' 




(K) 


(cm~^/km s~^) 


(cm-=») 



IN OUUO 


-L dust 


— 28 


< 11 25 




IC 860 


J- dust 


And 

— 4U 




0. i U It U.UZ 






onA -I- 7Q 

— zuu It /y 


±Z.44 It U. io 


4. / y It U.04 


NGC 5194 




— 40*^ 


< 10.67, 10.70 + 0.06 


< 4.89 


M 83 




= 31 


11.45 + 0.32 


5.29 + 0.10 




T/f (NH3) 


= 56+15 


11.49 + 0.40 


5.27 + 0.15 


Mrk 273 




= 48 


< 11.09, 11.13+0.12 


< 4.86 


NGC 5457 




= 40'' 


< 10.94 




IR 15107+0724 




= 40"^ 


12.46 + 0.12 


5.65 + 0.02 




Tif(NH3) -- 


= 189 ± 57 


12.71 + 0.12 


4.92 + 0.35 


Arp220 


Tdust 


= 44 


12.63 + 0.02 


5.64 + 0.02 




TjcCNHa) -. 


= 234 ± 52 


12.79 + 0.08 


4.09 + 0.09 


NGC 6240 




= 40'' 


< 10.61, 10.75 + 0.10 


< 4.80 


IR 17208-0014 


Tdust 


= 46 


< 11.05 




IR 17468+1320 




= 40'' 


< 11.11 




NGC 6701 


Tdust 


= 32 


< 10.63 




NGC 6921 


Tdust 


= 34 


< 10.65, 11.03 + 0.06 


< 4.27 


NGC 6946 


Tdust 


= 30 


11.01 + 0.36 


5.05 + 0.20 




Tk(NH3) 


= 47 + 8 


10.99 + 0.34 


5.06 + 0.21 


IC5179 




= 33 


< 12.28 




NGC 7331 




= 28 


< 10.81 




NGC 7479 


Tdust 


= 36 


< 11.03 




IR 23365+3604 




= 45 


< 10.89 




Mrk 331 




= 41 


< 11.05 





^ Tdz^st from Tablc[T] T/^(NH3) from |Mangum et al.|f2013[ l and |Ao et ar] | |201l[ | for NGC 1068. 
^ Sec g5.3| for assumptions used in calculating these quantities. 

^ LVG-dcrivcd column densities listed as limits for unconstrained density fits, followed by 
column density derived assuming a fixed density. See |5.3| for details. 
^ Assumed value. ^^^^ 
Kinetic temperature for the polar ring component (see |5.3.1| l. 




log(n(H„)) (cm-3) 



I 



Fig. 8. — LVG model predictions for the spatial density (n(H2)) 
and ortho-H2CO column density per velocity gradient (N(ortho- 
H2CO)/Av) for the 11 galaxies with a total of 14 velocity compo- 
nents in our sample with at least one detection of either the H2CO 
llO ~ 111 or 2ii — 2i2 transition where we have assumed that Tx = 
T^dust- Arrows indicate limits for galaxies with only one detected 
transition. Fit results are overlain on model line ratios (color) and 
intensities (solid (positive) and dashed (negative) contours, in K) 
for an assumed kinetic temperature of 40 K and no background 
continuum emission (Tc = 0). 



5.3.3. H2CO Excitation Pumping by Continuum Emission 

As was noted in [Mangum et~ar| ( |2008[ ), the H2CO 
lio ~ 111 transition can be driven into a maser excita- 
tion state by bright and compact nont hermal continuum 
emission at frequencies near 4.8 GHz ( Baan et al.||1986 ), 
producing anomalously bright H2CO emission. Analy- 
sis of the potential for maser amplification of the H2CO 
lio ^ 111 transition in Arp220, in light of the fact that 
the 2ii — 2i2 transition is observed in absorption, sug- 
gest either that the lio — In transition is not masing 
or that there are differences in the small-scale structure 
traced by the two K-doublet H2CO transitions. Using 
the simplest possible interpretation, we have so far relied 
on beam averaged quantities. While such an approach 
is reasonable for first order estimates, it ignores the po- 
tential impact of masers, distorting the picture outlined 
above. With our LVG models, all data could be suc- 
cessfully simulated without having to adopt population 
inversion. However, low-spatial resolution data alone do 
not provide a tool to discriminate safely between quasi- 
thermal and maser emission. Higher spatial resolution 
(i.e. using the VLA) measurements will make it possi- 
ble to clarify the influence of maser emission on the few 
extragalactic sources of H2CO K-doublct emission. 



5.4. H2CO Luminosity and Dense Gas Mass 

The correlation between the infrared continuum and 
molecular spectral line luminosity can be used to charac- 



HCN (Gao & Solomon 


2004a|b 


Carpio et ai. 


2UU6) extragalactic 



spatial density derived for our galaxy sample. 



star formation is the main power source for the large 
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Tj = T,t(NH3) 



Arp 220 IR 15107+0724 

(Tg = 234±52 K) (Tj = 189±57 K) 




LVG Model-Derived 

Density and HgCO Column Density 



NGC 3079C1 

(Tj > 100 K) 



NGC 3079C2 

(Tj > 150 K) 



M 83 

'(\ = 56±15 K) 



NGC 660 

(Tj = = 37 K) 



6 



log(n(H2)) (cm-^') 



Fig. 9. — LVG model predictions for the spatial density (n(H2)) and ortho-H2CO column density per velocity gradient (N(ortho- 
H2CO)/Av) for the 12 galaxies and galaxy velocity components in our sample with measured H2CO 1 10 — 111 and 2ii — 2i2 and NHs-derived 
Tfc (Mangum ct al. 2013). For each galaxy or galaxy component the NHs-measured T/^ used in the LVG model fit is listed. Galaxies with 
a lower limit to the kinetic temperature are displayed as a range in (n.N) at the lower-limit kinetic temperature, with an arrow indicating 
the trend in (n,N) as the kinetic temperature is raised. As the fit range is so large for Maffei2Cl (due to low signal-to-noise) and the 
high-temperature fit to N GC 66 0, the fit ranges for these galaxy components are shown dotted. The need for the two fit ranges shown for 
NGC 660 is described in fjOTT] 



infrared luminosities observed in these galaxies. Fol 
lowing the formalism for calculating molecular spectral 
line lumino sity and dense gas mass described in Mangum 
et al. I ( 12008 ), we have derived Lir ( [Mangum et al.||2008 
Equation 12) fr om existing infrared measurements anc 
Mdense(H2C0) ( [Maugum et al. 2008, Equation 11) from 
our H2CO measurements (I'ablc [7]) : 



M dense — ^''mol^'-s^A 



1.125 X 10*^ TT0'^^{arcsec)x 

Dl{Mpc)N,noiicm-^)^mH2 jgm) 
41n(2)(l + z)2X™o, 
5 . 79 6*2 {arcsec) Dl (Mpc) N^^i ( x 10i°cm- ^ ^ 



(l + z)2X„„i(xlO-9) 



Ma 



(2) 



with the gaussian telescope beam, Dj^ is the angular size 
distance to the galaxy, to//^ is the mass of molecular 
hydrogen, fi accounts for the mass fraction due to He 
(1.36), X,noi is the abundance (relative to H2) of the 
molecule, and uniform filling of the molecular volume is 
assumed. Where two calculations of Mdense(H2C0) are 
listed the two calculations of N(ortho-H2CO)/Av listed 
in Table |6] have been used. 

We also list M^ensefHCN) derived from HCN mea- 
surements (Gao & Solomon 2004a) . Note that our 
Mdense(H2C0) Calculations assume X(H2C0) = 10^^ 
and include a conservative 50% uncertainty to our as- 
sumed source sizes and H2CO abundances. Due to the 
lack of H2CO source structure information our H2CO- 
derived Mdense values are only estimates to be compared 
with other molecule-derived dense gas masses. 

Our Mdense estimates range from < 4 x IO^Mq 
(NGC 598) to - 1 X lOi°M0 (NGC 660, NGC 1144, 
IR15107-f0724, NGC 6921). In general, Mdense derived 



where fls is the solid angle of a gaussian source convolved from our H2CO measurements agrees to within an order- 
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Infrared Luminosities and Dense Gas Masses 



Galaxy 


Sizc^ 


T b 


Mde„ae(HCN) = 
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^ Sec j ]5.1| for galaxies detected in H2CO. Otlier source size reference s and measu rements used 
arc NGC55, NGC 1022 NGC 1055, NGC 3627 , NGC 5005, NGC 6701: |Gao fc Solo mon |20 04all, 
HCN J=l - 0; IC1623: [imanishi et al.| l |2007] , H CN and HCO+ J=l - 0; NGC 598: |Rosolowsky| 
|et aLlpOOT} , CO J=l - ; NGC i?2S7NUUl736: |Leroy et al.| l |2009| l, CO J=2 - 1; LJGC0;!369: 
this work; IR 03359+1523: 'Sanders, Scoyille, fc Soifer| l |199'T^ 7CO J = l - 0; NGC 1614, NGC 2623, 
Arp 55, IR 17208-001 4: lono et al. (2009|, CO J=3 - 2 and HCO+ J=4 - 3; yiIZw 31, NGC 3690: 
[Gracia-Car pio et al.* ('20081, HCN ,1 = 1 - 0; IR 10173+0828: iPlanesas et al.|l|1991( , CO J=l - 0; 
NnrrS^TTT Liscnfcld et al. (2008), CO ,1=1 - 0; NGC 4631: |lsraelHiiOOiH l, Cc3 J =2- 1 and 3-2; 
NGC 5457: iHelfer et al.,(,2003J, CO J=l - 0. 

^ Luminosities from |Sande^et al. 1 | |2003i, derived frorn IRAS fluxes over 8 to 1000 am. 

Except for Arp 220 and NGC 6240 ) |Greve et al.|2009l , from [Gao &: SolomonH2004a| |. 
^ Limits assume 3cr in nondetected H2CO transition(s). 

^ When two values for M^e^ge ( H2CO ) are listed they are derived from the two calculations of 
N(ortho-H2CO)/Av described in ^5.3| for galaxies with only one H2CO transition. 
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Fig. 10.— LVG model H2CO lio - In and 2ii - 2i2 line ratio 
(color) and H2CO lio — In intensity (solid (positive) and dashed 
(negative) contours, in K) predictions as functions of the spatial 
density (n(H2)) and ortho-H2CO column density per velocity gra- 
dient (N(ortho-H2CO)/Av). LVG model predictions at kinetic 
temperatures of 50 K (top) and 200 K (bottom) are shown (Tc 
= 0). Note how the predicted spatial density for a given H2CO 
line ratio decreases (increases) as the kinetic temperature increases 
(decreases). 

of-magnitude with Mdense derived from HCN measure- 
ments (excluding NGC1144, which differs by a factor 
of ~ 100). There is an apparent trend, though, for 
Mdense(H2C0) < M^ense (HCN) . Notc that the uncer- 
tainties in many of the quantities which go into the cal- 
culation of Mdense limit its accuracy to no better than 
a n order-o f-magnitu de^ 

found that the values of 



10" 



NGC214a^ I 



10" 



10° 



Fig. 11. — L[ji versus M^^ng,, for the 11 galaxies and one galaxy 
offset position (M 82 and M 82SW overlap) where we have derived 
values for N(ortho-H2CO)/Av from our LVG model fits, with 
derived from NH3 measurements (Mangum ct al. 2013|l. M^je^ge 
has been derived using Equation [2] assuming X(H2CO) = 10~^ 
and uncertainties calculated assummg 50% errors in source size and 
X(H2CO). Filled circles and squares are used to indicate H2CO- 
detected galaxies at distances > and < 50 Mpc, respectively. Limit 
arrows indicate lower-limits which extend to zero. The dotted line 
represents the log-linear relation \og{Ljji) = 2 + 1.0 log(A/^i3„se) 
and is not a fit to these data. 



Mdense (HCN) are factors of 5-200 smaller than M{H2) 
derived from CO measurements for their sample of spi- 
ral, LIRG, and ULIRG galaxies, many of which are in- 
cluded in our sample. Our Mjjense (H2CO) values are 
generally smaller than those derived using HCN. The 
progression from high to low masses derived from dense 
gas tracers such as CO, HCN, and H2CO reflects the 
hierarchical str ucture of the giant molecular clouds in 
these galaxies (Gao & Solomon 2004a|b Greve et al. 
2009 1. H2CO traces a denser, more compact, component 
of the giant molecular clouds in our galaxy sample than 
low-excitation transitions of CO or HCN. This result is 
consistent with high-resolution stu dies of the K- double t 
H2CO emiss ion in our own Galaxy (Mangum et al.|1993[ ) 



In Figure 11 



wc show Mdense versus Lm for the 11 
galaxies and one galaxy offset position where we have 
derived values for N(ortho-H2CO)/Av from our LVG 
model fits (rather than simply limits to the ortho-H2CO 
column density) with Tj^ d erived from NH3 



measure- 



11 



ments ( Mangum et al.[2013 ). We also show in Figure 
the log-linear relation log(i7i^) = 2-1- l.Q\og{Mdent 
Note that given the uncertainties regarding the size of 
the dense gas distribution wc have not performed a for- 
mal fit of this relation. This correl ation, first noted by 



Solomon, Downes, & Radford ( 1992 ), reflects the correla- 
tion between the infrared luminosity and the amount of 
material to form stars in galaxies. This in turn leads 
directly to the Kennicutt-Schmidt laws, which relate 
the star formation rate to the mass of gas available to 
produce stars (Schmidt 1959] |1963 1 in galaxies. Since 



Gao & Solomon (2004a 



Mdense oc Lmol , thc suggcstiou ot a linear correlation be 
tween Ljji and Mdense{H2C0) implies that, similar to 
HCN, H2CO traces the dense star- forming gas in star- 
burst galaxies. 
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As was alluded to in Mangum et al. ( 2008| p the dense 
gas mass can also be calculated using the spatial density 
n(H2) rather than the molecular column density l^moi- 



dense 



n{H2)hQ.sD\m,-aoi 



3.47 X 10^^ Tieli arcsec) x 

Dl{Mpc)n(H2){cm^^)h{pc)iJ,mH2{9'm) 
41n(2)(l + z)2 

1.31 el{arcsec)D'l{Mpc)n{H2){cm-'^)h(pc) 

(1 + Z)2 



Mo, 



Mo 



(3) 



where we have assumed that the volume of emitting gas 
is a uniformly-filled face-on disk with height fi = 1.36, 
and rrimoi is the mass of an individual molecule. Com- 
paring the mass derived from Equation [21 (Mdense) with 
M' 





TABLE 8 




Comparison of Mdense and M'^^^ 




Galaxy 


dense {"-2^^) 


Mj CHoCO)'' 

dens e \ ^ j 




(10 Mq) 


do* A/r„ "I 


NGC253 {Tk = 78±22K) 


9.68±6.97 


12.68±6.57 


NGC660 (Tk = 160±97K) 


82.04±124.62 


45.12±123.47 


Maffoi 2 


0.33±0.25 


3.31±2.18 


IC342 (T/f = 24 ± 7K) 


0.16±0.16 


4.05±3.37 


NGC2146 


6.19±4.70 


105.22±55.59 


M82 


19.10±13.76 


31.21±17.47 


M82SW 


18.68±13.66 


42.10±26.53 


NGC 3628 


2.78±2.07 


12.90±10.71 


M83 


1.77±2.25 


36.39±31.65 


IR 15107+0724 


110.00±83.65 


49.92±46.51 


Arp220 


62.85±45.94 


1.46±0.96 


NGC 6946 


0.12±0.13 


3.11±3.50 



Derived from Equation 
Derived from EquationjS] 



dense 



allows for a consistency check oY our dense gas emission or absorption lies within a rather limited range 



mass calculation. Comparing the terms in Equations [2] 
and [3j this consistency check amounts to a check of our 
assumptions regarding the H2CO abundance (X„jo/) and 
the dense gas emitting volume thickness h. Assuming 
ft. = 10 pc (similar in size to the CMC structures iden- 
tified in IC342 by [Meier et al.| ( [201 1[ )) for all galaxies 
we find that for the ten galaxies winch have both N,„oi 
and n(H2) measurements that our calculations of M^g^^^ 
and Mdense (Table [sl agree to within a factor of ~ 2 for 
NGC 253, NGC6607^M82, and IR 15107+0724. For the 
other seven galaxies MJ^^ and Mdense agree to within 
a factor of ~ 20, with Md^^^^ larger than Mdense in aU 
but one galaxy (Arp220). To bring these two dense 
gas mass estimates into agreement we can either lower 
the dense gas emitting volume thickness h or lower the 
H2CO abundance by a factor of ~ 20. As our estimate 
for h = 10 pc is consistent with the GMC-scale struc- 
tures imaged at high-spatial re solution in several nearby 
starburst galaxies (e.g. IC 342 ; [MeieretalJ M82; 



of 10 



4.5 



Carlstrom & Kronberg ( 1991 )), the more plausible option 
appears to be to lower the H2CO abundance to ^ 10~^° 
(recall that our assumed X,„o/ for all calculations above 
is 10~^). Note that models of molecular abundances in 
starburst galaxies with low to moderate cosmic ray ion- 
ization rates (the local Milky Way value of 2 x 10~^^ 

to 5 X 10-^6 ! 
10-10 _ 10-9 

jerink et al.[ [2' 



) predict H2CO abundances in the range 
jBayet et al.|201l[) and 10"^ - lO^^^ (jMei 
Oll[ : R. Meiierinfc 



Meijerink 2011, priv. comm.). In 
these models as the cosmic ray ionization rate is increased 
the H2CO abundance decreases . If mechanical heating is 
included ( Meijerink et al.|2011 ) a more complicated pre- 
diction 01 the H2CU abundance with increasing cosmic 
ray ionization rate is predicted which lies in the range 

io-Mo-13. 

6. DENSITY-INDEPENDENT STAR FORMATION 



As was noted in ^ 5.3.2 



the measured mean spatial den- 
sity within the starburst environments of our sample of 
13 galaxies with measured H2CO lio — lii and 2ii — 2i2 

® Equation 14 in |Mangum et ar] | |2008| presents a similar relation 
for M^^^^^ that has two errors: a missing factor of ^i^^yi an 

incorrect exponent to D^. Equation 14 in |Mangum et aL] ( |2008| is 
also not properly normalized. 



10 cm These galaxies span a consider- 



able variety in kinetic temperature (Tk — 30 — 300 K) 
infrared luminosity {Lju ~ 10^'^ — and dense 

gas mass {Mdense — 10^ — 10^^ Mq), representing a wide 
range of starburst environments. The physical size scales 
over which these measurements apply range from < 0.3 
to ^ 20 kpc based on our spatial resolution of 51" to 153", 
sample galaxy distances of 1 to 200 Mpc (see Table [I]) 
and unknown single dish beam filling factors. The rel- 
atively narrow range of the mean spatial density within 
our starburst galaxy sample suggests that spatial density 
is not a driver of the star formation process in the most 
luminous starburst galaxies. Furthermore, a relatively 
constant mean spatial density implies that the Schmidt- 
Kennicutt relation between L/^ and Mdense is a measure 
of the dense gas mass reservoir available to form stars, 
and not a refiection of a higher spatial density in the 
most luminous and massive starburst galaxies. 

We should point out that the relatively narrow range 
of spatial densities derived in our starburst galaxy sam- 
ple is not due to limitations of our densitometry tech- 
nique. The H2CO K-doublet transitions we have used in 
this analysis absorb the cosmic microwave background 
at n(iJ2) ^ lO'^ '^ cm~'^ and appear in emission at higher 
spatial densities (see For typical II2CO abundances 
the II2CO K-doublet transitions are excited at n{H2) 
> 10^ cm~'^. As there is no upper limit to the spatial 
density sensitivity of these H2CO transitions, biased sen- 
sitivity to spatial density does not affect our spatial den- 
sity measurements. Observational support for this lack 
of bias to spatial density e xists for both lower dZeiger fc 

^ - - ' inrr — 



201 



Gms- 



Darling 2010) and higher (McCauley et al 
Eurg et al.||201ip spatial densities than those derived in 
this study. 

7. CONCLUSIONS 

Measurements of the lio — In and 2ii — 2i2 K-doublct 
transitions of II2CO toward a sample of 56 starburst 
galaxies have been used to study the physical conditions 
within the active starburst environments of these galax- 
ies. By applying detections of at least one of these K- 
doublet II2CO transitions toward 21 of these starburst 
galaxies an estimate to the spatial density (n(H2)) and 
dense gas mass within each of these starburst galaxies has 
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Mangum ct al. 



been derived. For the 13 galaxies where both K-doublet 
H2CO transitions have been detected we have derived 
accurate measurements of the mean spatial density. By 
applying an appropriate measurement of the dense gas 
kinetic temperature in 11 of the galaxies in our sample 
we have improved t he spatial densi ty measurement pro- 
cedure presented in Mangum et al. (2008). Furthermore, 
all of our K-doublet li2CJU measurements can be fit to 
our LVG model assuming quasithermal exciation. This 
fact does not enti rely exclude po tential maser emission, 
as noted in ^Mangum et al. ( 2008 1 . 

Using measured kinetic temperatures has narrowed the 
range of de rived spatial densities presented in IManguml 



et al. (2008 ) to lO''-^ to 10^-^ cm ^ in our starburst galaxy 



sample! This improvement to the spatial density mea- 
surement in our starburst galaxy sample has removed 
the trend between Lj j^ and our derived n{H2) noted in 
Mangum et al. (2008). That trend, which suggested that 
starburst galaxies with higher spatial densities also pos- 
sessed higher IR luminosities, was indicative of another 
representation of the Lm-Mdense correlation. Instead, 
our results now imply that the Schmidt-Kennicutt rela- 
tion between Ijir and yifi^nse'- 

• Is a measure of the dense gas mass reservoir avail- 
able to form stars, and 

• Is not directly dependent upon a higher average 
density driving the star formation process in the 
most luminous starburst galaxies. 

This extension of the characterization of the spatial den- 
sity in starburst galaxies presented in |Mangum et al. 
p008J produces the most accurate measurements of this 
important physical quantity in starburst galaxies made 
to-date. 



As was done in Mangum et al. (2008), we have also 



used our H2CO measurements to derive a measure of the 
dense gas mass which ranges from 0.14— llOx lO^M©. By 
comparing this traditional measure of the dense gas mass 
to that derived using our spatial density measurements 
we find agreement to within a factor of ^ 20. The most 
extreme differences between these two methods of calcu- 
lating the dense gas mass using H2CO can best be rec- 
onciled with a modification of the assumed H2CO abun- 
dance of X(H2C0) = 10"^ toward several of the galaxies 
in our sample. Recent modeling of the molecular abun- 
dances in starburst galaxies sugg est that low to moderate 
cosmic ray ioniza tion ra tes (Bay et et al. 2011) and me- 
chanical heating (Meijerink et aE' 2011 1 can affect the 
H2CO abundance in starburst environments. Further- 
more, comparison of our measures of the dense gas mass 
using H2CO with those derived using HCN suggest that 
the H2CO K-doublet transitions trace a denser, more 
compact, component of the giant molecular clouds in our 
starburst galaxy sample than low-excitation transitions 
of CO or HCN. 
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